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FOREWORD 


This  study  was  initiated  by  the  Reconnaissance  Charting  Branch  of 
the  Intelligence  Laboratory  as  a  part  of  the  Air  Force  program  to  develop 
a  target  location  and  charting  capability  compatible  with  future  weapon  systems. 
Specifically,  the  purpose  of  this  study  was  to  investigate  the  unique  possibilities 
presented  by  the  inclusion  of  different  kinds  of  sensors  in  a  high-performance 
boost-glide  type  of  vehicle.  However,  to  keep  the  work  within  the  limitations 
of  the  contract,  only  the  various  combinations  of  panoramic  photography  and 
side  scan  radar  were  considered.  Because  of  the  lack  of  the  fundamental 
experiments  that  are  needed  for  an  adequate  determination  of  the  expected 
atmospheric  effects,  that  portion  of  the  work  was  only  briefly  covered. 

Two  important  conclusions  were  drawn  from  this  work: 

1,  To  match  different  sensor  records  from  the  high  altitudes 
under  consideration,  the  establishment  of  a  stereo  model  or  even  the  rectification 
of  the  images  would  not  permit  registration  of  the  images,  because  of  the  dif¬ 
ference  in  direction  of  relief  displacement,  unless  a  local  datum  and  techniques 
for  forced  matching  of  recognizable  points  were  used. 

2,  Before  a  choice  is  made  for  a  multi-sensor  collection 
system,  many  factors  need  to  be  considered,  e.  g,  ,  since  radar  is  essentially 
independent  of  stabilization  and  radar  stereo  is  only  slightly  dependent  upon  the 
relative  positions  of  the  two  sensors,  a  radar  system  will  outperform  the  photo¬ 
graphic  system  if  navigation  and  stabilization  errors  are  large  whereas,  at  low 
altitudes,  stabilization  errors  do  not  appreciably  affect  the  systems  and  photo¬ 
graphic  sensors  can  be  used.  At  higher  altitudes,  stabilization  errors  can 
dominate  a  photographic  system.  Implied  in  these  conclusions  is  the  fact  that 
resolution  is  not  a  limiting  factor. 

A  few  conditions  for  each  parametric  equation  were  used  as  a  basis 
for  the  above  conclusions  as  presented  in  this  report.  Because  it  was  not  advisable 
to  examine  all  the  possible  conditions,  those  of  greatest  interest  to  an  individual 
may  not  have  been  included.  If  the  necessary  information  is  not  found  in  this 
document,  it  may  be  possible  to  utilize  the  computer  program  to  generate  the 
information. 
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This  work  has  resulted  in  interesting  conclusions  that  should  be 
taken  into  consideration  during  the  configuration  of  a  vehicle  to  be  used  for 
multi- sensor  target  location  or  charting. 

This  document  presents  the  results  of  a  scientific  investigation.  It 
is  hoped  that  it  will  stimulate  interest  in  the  problem  under  consideration  and 
will  serve  as  an  indication  of  our  interest.  Comments  on  this  work  are 
invited  and  should  be  addressed  to  the  contract  scientist. 


Lynwood  D.  Sinnamon,  Jr. 

Physical  Scientist 

ATTN:  RCWICC 

Rome  Air  Development  Center 

Griffiss  AFB,  New  York 

March  1961 
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ABSTRACT 


Target  location  and  charting  using  data  acquired  by  boost-glide 
type  vehicles  are  investigated.  Particular  emphasis  is  placed  on  two  aspects 
of  the  problem;  the  use  of  multiple  sensors,  and  atmospheric  effects. 

The  possibility  of  Improving  target  location  accuracy  by  using 
combinations  of  different  sensor  types  is  investigated.  A  parametric  analysis 
provides  error  contributions  due  to  various  sources  as  functions  of  altitude, 
target  offset,  and  sensor  separation  in  3-D  location  systems.  Trade-offs 
between  various  systems  are  presented. 

Atmospheric  effects  are  described  analytically.  Where  experi¬ 
mental  data  could  be  obtained,  various  refraction  effects  are  evaluated 
numerically  and  plotted  as  a  function  of  altitude. 
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SUMMARY 


A  one  year  study  has  been  conducted  during  which  two  aspects  of 
the  problem  of  target  location  and  charting  from  boost-glide  type  vehicles  have 
been  investigated.  First,  possible  improvements  in  target  Location  accuracy 
through  the  use  of  combinations  of  different  sensor  types  were  studied.  Second, 
the  atmosphere,  a  major  source  of  image  distortion  and  degradation,  has  been 
treated  analytically,  and  a  survey  of  present  experimental  data  was  conducted. 

Target  location  accuracy  was  investigated  by  evaluating  location 
errors  for  several  combinations  of  sensors.  Sensor  altitude,  target  position, 
and  sensor  separation  were  varied  parametrically.  Measurement  errors  were 
varied,  so  that  system  trade-offs  could  be  found  as  a  function  of  improvement 
in  sensors  or  other  equipment  (such  as  navigation  systems).  It  is  likely  that  a 
combination  of  radar  and  camera  in  a  single  vehicle  can  outperform  either 
radar  or  photographic  stereo  target  location  systems  singly,  unless  the  base 
line  for  photographic  triangulation  is  large  (as  in  the  case  of  convergent  photo¬ 
graphy  or  photographs  taken  on  separate  flights,  where  side-lap  is  used). 

Radar  systems  are  less  affected  by  navigation  and  orientation  errors  than  canriera 
systems,  and  so  improve  relative  to  photographic  methods  as  these  errors 
increase. 

As  a  by-product  of  this  study,  it  was  possible  to  estimate  the  accuracy 
with  which  images  could  be  brought  into  registration  for  superposition.  It  was 
found  that  single  sensor  images  contain  relief  displacements  large  enough  to 
allow  only  coarse  registration.  Matching  stereo  models  improves  the  situation, 
but  auxiliary  aids  such  as  relative  control  points  must  be  used  if  small  target 
objects  are  to  be  registered. 

Atmospheric  refraction  studies  indicate  that  above  100,000  feet, 
in  the  boost-glide  altitude  regime,  the  major  distortion  and  degradation  is 
due  to  atmospheric  turbulence.  Present  experimental  data  are  insufficient  to 
establish  reliable  values  for  parameters  in  the  theoretical  equations,  so  precise 
evaluation  of  atmospheric  effects  is  not  possible  at  this  time. 
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It  is  recommended  that  further  work  be  done  on  target  location 
•yateme  to  determine  the  reduction  of  errors  when  least  squares  adjustment 
of  data  is  used  in  conjunction  with  ground  control.  Atmospheric  experi¬ 
mentation,  especially  in  the  area  of  turbulence,  should  be  continued. 


1.  INTRODUCTION 


The  Cornell  Aeronautical  Laboratory,  Inc.  ,  under  contract 
AF30  (602)-2183  with  the  Rome  Air  Development  Center,  has  studied  problems 
of  reconnaissance  from  boost-glide  type  vehicles.  The  contract  covered  one 
year  at  approximately  a  one  man  level  of  effort. 

Such  a  study  has  many  facets.  Two  choices  are  open  if  a 
relatively  small  study  Is  to  be  conducted.  Either  many  aspects  may  be 
investigated  in  cursory  fashion,  or  a  few  may  be  more  thoroughly  studied. 

With  the  agreement  of  the  sponsoring  agency  this  Laboratory  chose  the  second 
alternative,  with  study  effort  being  concentrated  on  the  problems  of  target  location 
tion  using  multiple  sensors  and  atmospheric  optical  effects. 

The  question  which  arises  with  respect  to  target  location  and  charting 
is:  Can  point  location  accuracy  be  improved  by  using  combinations  of  sensors 
for  location  purposes?  It  is  reasonable  to  expect  that  a  combination  of  camera 
and  radar  might  outperform  either  sensor  used  separately,  since  the  camera 
provides  directional  data  while  the  radar  yields  accurate  range  measurements. 

The  target  location  accuracy  study  described  in  this  report  shows  the  conditions 
under  which  each  sensor  combination  can  be  used  to  best  advantage.  The  problems 
of  mensuration  are  considered,  with  the  assumption  that  targets  are  identifiable. 
Although  the  identification  problem  is  not  directly  treated,  this  same  study  also 
provides  data  concerning  the  possibility  of  superimposing  images  from  different 
sensors. 

The  atmosphere  is,  of  course,  a  major  source  of  image  degradation 
and  distortion  in  boost-glide  reconnaissance.  The  study  of  the  atmosphere  is 
essentially  a  review  of  the  state-of-the-art  of  this  field.  Refraction  is  treated 
in  terms  of  the  various  sources  of  steady  and  fluctuating  density  gradients  in 
the  atmosphere.  Results  are  j>resented,  using  the  best  available  experimental 
data,  which  indicate  the  relative  effect  of  these  sources  of  image  degradation 
and  distortion. 

The  body  of  the  report  is  divided  into  two  parts.  Part  1  describes 
the  target  location  study  while  Part  11  is  concerned  with  the  atmosphere. 
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PART  1 


2.  TARGET  LOCATION 


Target  location  may  be  considered  in  several  stages  of  sophistication, 
The  simplest  form  of  target  location  consists  of  using  a  single  sensor  to  obtain 
an  image  which  contains  the  target.  The  range,  or  the  direction,  to  the  target 
is  given  by  measurements  of  this  image.  The  sensor  altitude  is  then  combined 
with  these  measurements  to  compute  the  target  location  with  respect  to  the 
sensor.  A  knowledge  of  the  sensor  position  with  respect  to  the  base  of  operations 
completes  the  target  location.  This  simple  system  is  plagued  by  two  major 
errors,  both  due  to  the  earth  reference  used.  The  simplest  calculations  assume 
a  flat  earth,  and  there  is  no  way  to  account  for  terrain  altitude  variations  in  any 
ease.  The  effects  of  earth  curvature  and  terrain  relief  are  calculated  in  the 
following  section.  Terrain  relief  errors  alone  are  sufficiently  large  to 
warrant  study  of  means  to  minimize  their  effect. 

The  errors  due  to  using  a  reference  datum  can  be  eliminated  by 
using  data  from  two  sensor  images  rather  than  one.  "Stereo"  target  location 
relies  only  on  measurements  on  sensor  images  and  a  knowledge  of  the  sensor 
position  with  respect  to  the  base.  Each  sensor  provides  one  or  two  measurements 
which  aid  in  locating  the  target.  Therefore,  since  the  target  must  be  located 
with  respect  to  three  coordinates »  two  sensor  Images  are  needed.  The  mathe- 
matics  of  target  location  by  this  means,  with  an  error  analysis,  is  the  subject 
of  the  major  portion  of  Part  1  of  this  report.  Various  combinations  of  sensors 
are  considered,  indicating  which  combination  provides  the  most  accurate 
target  location  under  various  circumstances. 


"Stereo,  "  as  used  in  this  report,  means  any  system  which  locates  targets  by 
making  use  of  two  sensor  records  (at  least)  to  provide  three-dimensional  loca¬ 
tion  information.  The  term  is  not  restricted  to  systems  in  which  the  Images 
are  viewed  as  a  stereo  model.  Although  this  is  3-D  target  location,  the  term 
"stereo"  was  chosen  to  eliminate  confusion  between  "3-D  radar"  which  uses  one 
sensor  and  "Stereo  radar"  which  employs  trilate ration. 

ee 

One  exception  to  this  is  3-D  radar,  in  which  range,  azimuth  and  elevation  are 
all  measured.  However,  elevation  measurement  errors  are  too  large  to 
seriously  consider  this  system  for  target  location  from  high  altitude. 
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Any  measurementB  which  are  available  in  excess  of  the  required 
three  may  be  used  to  reduce  sensor  measurement  errors  by  least  squares 
fitting  of  data.  Thus,  if  one  target  appears  on  several  images  (such  as 
happens  in  panoramic  photography  at  large  offset)  the  error  expected  from 
a  simpler  analysis  could  be  reduced  somewhat.  Navigation  errors  can  be 
reduced  if  several  points  on  successive  images  can  be  used  to  extend 
ground  control.  Relative  control  points  would  reduce  error  contributions 
due  to  navigation  measurements  even  if  no  absolute  ground  control  was  avail¬ 
able.  Again,  this  error  reduction  would  involve  a  least  squares  fitting  procedure. 
Such  extensive  mathematical  analyses  are  beyond  the  scope  of  the  present  report. 
However,  an  introduction  to  this  subject,  showing  that  such  procedures  as  have 
been  developed  for  photography  can  be  extended  to  include  other  types  of  sensors, 
is  presented  as  an  appendix  to  this  report. 
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3.  ERROR  PROPAGATION  IN  NONSTEREO  SYSTEMS 


Two  major  errors  arise  when  an  attempt  is  made  to  locate  a  target  by 
using  a  single  sensor  and  a  reference  datum  surface.  The  first  contribution, 
which  is  negligible  at  low  altitude,  is  the  effect  of  earth  curvature.  This  may, 
of  course,  be  eliminated  by  the  use  of  a  curved  reference  datum.  This  effect 
is  considered  here  only  to  the  extent  required  to  indicate  when  the  simpler 
flat  datum  may  be  used. 

The  second  error,  similar  in  form  to  the  first,  is  the  effect  of  terrain 
relief.  The  location  error  due  to  terrain  relief  will  be  evaluated. 


3.  1  Earth  Curvature  Effect 

Figure  3.  1  illustrates  the  geometry  of  the  situation  for  both  camera 
and  radar.  The  error  in  using  a  flat  earth  assumption  can  be  computed  as: 

Error  »  Fa  »  NP~NP^  •  camera  3-1 

Efi  “  ,  for  radar  3-2 

The  distances  HP,  NP\  HP'*  can  be  computed  from  the  figure,  as: 

3-3 


NP  -CC /St 


eC  ■ 


For  the  radar  case: 


Since  ^  is  an  angle  of  a  few  degrees  or  less,  3-6  may  be  written: 
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sensor 


-  radius  of  earth 
H  s  sensor  altitude 

N  =  nadir  point 

P  s  target  point 

P*  =  apparent  target  location  (camera,  flat  datum) 

s  apparent  target  location  (radar,  flat  daitum) 

^  =  field  angle 


Figure  3.  1  EFFECT  OF  DATUM  SURFACE  ON  TARGET  LOCATION 
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Figure  3.  2  shows  the  error  in  assuming  a  flat  earth  datum  as  a 
function  of  angle  of  the  target  away  from  nadir.  Several  altitudes  are  shown. 
Only  errors  for  the  photographic  system  are  plotted,  since  Equation  3-7 
shows  that  the  error  in  the  radar  case  is  essentially  a  linear  function  of 
ground  offset.  It  should  be  noted  that  photographic  and  radar  effects  are 
opposite  in  sign. 


3.  2  Terradn  Relief 

Figure  3.  3  shows  the  effect  of  terrain  relief.  The  pertinent 
equations  relating  to  target  location  are  given  below.  (It  is  assumed  that 
the  target  is  at  a  position  at  right  angles  to  the  flight  line,  directly  to  the 
side  of  the  sensor  platform. )  The  Y  position  of  the  target  is  computed  as: 


for  radar 

3-8 

Yt-- 

for  camera 

3-9 

Relating  these  quantities,  using  Figure  3.3,  assuming  the  target  to  be  a 
height  h  above  the  terrain  reference  under  the  sensor: 

*  '■ISI  I  II  • 


s  * 


(/✓-a; 


4 


(H-h)'- 


-/ 


3-10 


sensor  position 


FIGURE  3.  3  Effect  of  Terrain  Rf-llef  on  Target  Location 
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A  numerical  example  can  be  presented  to  indicate  possible  values 
of  terrain  relief  errors.  An  angle  within  the  field  of  both  the  camera  and  radar 
was  chosen.  The  relief  displacement  for  radar  images  increases  for  smaller 
angles  away  from  nadir,  while  for  photography  the  relief  displacement  is 
greater  at  larger  angles.  Table  3.  1  summarizes  the  results. 

Table  3.  1  RELIEF  DISPLACEMENT  ERROR 
(Using  da^um  reference  beneath  sensor) 


sensor  altitude 

angle  from  nadir 

relief 

actual  offset 

relief  displacement 

H 

Photo 

Radar 

100  Kft 

25* 

1  Kft 

46  Kft 

+  466  ft 

-2140  ft 

100  Kft 

25* 

100  Kft 

46  Kft 

+  46.  6  ft 

-  214  ft 

100  Kft 

6* 

1  Kft 

10  Kft 

+  105  ft 

(not 

100  Kft 

6* 

100  Kft 

10  Kft 

+  10.  5  ft 

applicable) 

3.  3  Results 

It  can  be  seen  from  Sections  3.  1  and  3.  Z  that  accurate  target 
location  cannot  be  achieved  with  nonstereo  systems.  Since  the  use  of  stereo 
would  reduce  target  location  errors,  stereo  systems  will  be  considered  in 
following  sections. 

If  an  attempt  were  made  to  bring  a  complete  radar  image  and  a  photo¬ 
graph  into  registration,  the  images  of  some  points  would  match,  while  others 
would  be  displaced  relative  to  each  otlier  by  the  order  of  thousands  of  feet. 
Clearly  only  very  large  targets  would  still  register  under  these  conditions. 
Since  large  targets  are  already  usually  detectable  on  both  images  separately, 
there  would  be  little  advantage  in  combining  them.  Improvement  could  be 
achieved  in  regions  where  identifiable  points  allowed  a  local  datum  and  scale 
to  be  used. 

Assuming  that  only  narrow  field  photographs  are  to  be  registered, 
the  situation  is  better.  If  a  target  is  near  the  center  of  one  photograph,  no 
appreciable  relief  displacement  occurs.  If  the  same  target  appears  on  a 
second  photograph  at  6*  off  axis.  Table  3.  1  shows  that  the  target  may  be 
shifted  ten  to  a  hundred  feet  due  to  relief.  If  reliefs  are  in  the  order  of  100  ft 
(gently  rolling  terrain),  the  10  ft  displacement  is  typical,  in  which  case  some 
value  may  he  derived  from  superposition. 
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4.  STEREO  TARGET  LOCATION  ERROR  ANALYSIS 


Aa  a  basis  for  the  parametric  study  to  be  considered,  this  section 
describes  the  basic  target  location  equations  for  various  sensor  combinations 
and  the  error  analysis  which  uses  these  equations.  The  results  of  this  section 
are  used  in  the  following  section  which  actually  describes  the  parametric  study. 

The  equations  will  be  developed,  but  not  applied,  in  the  present  section. 

It  is  assumed  that  a  single  target  is  to  be  located  with  respect  to  a  coordinate 
system  fixed  at  a  known  base.  The  target  is  to  be  located  in  all  three  dimensions, 
that  is  latitude,  longitude  and  altitude  (or  equivalently  t, /,  e  in  a  Cartesian 
system).  Since  no  single  sensor  in  one  position  provides  the  necessary  data 
for  3-0  target  location,  the  problem  must  involve  two  sensors.  These  might  be 
either  the  same  sensor  at  different  locations,  or  two  different  sensors.  The 
problem,  then  is  to  determine  a  vector  from  the  coordinate  origin  at  a  base 
to  the  target,  by  using  partial  information  about  two  vectors  from  sensors  to 
the  target  plus  navigation  vectors  from  the  base  to  the  sensors  and  orientation 
information  about  the  sensors.  Figure  4.  1  illustrates  the  problem.  The  vector, 

P  ,  is  to  be  determined. 

Generally,  P  is  determined  by  solving  the  following  vector  equation; 

P^S,*K,  4-1 

where  the  navigation  vector,  5^  is  known  in  base  coordinates  and  must 
be  solved  for  by  using: 

R,  ^  4-2 

and  ?  -  5^  -t,  4-3 

In  equation  (4-2),  only  partial  information  about  ^  and  is  available 
from  measurements.  Therefore,  the  equations  implied  by  (4-2)  actually 
differ  in  form  for  each  combination  of  sensors  considered,  and  will  be  set 
up  and  solved  separately  for  each  case.  If  the  two  sensor  stations  are  part 
of  a  single  flight  over  the  area,  S,  the  relative  position  of  station  2  with 
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base  system 


Figure  4,  1 


STEREO  TARGET  LOCATION 
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respect  to  1,  f  s  in  error  only  due  to  navigation  from  station  1  to  2.  On  the 
other  hand,  if  the  two  sets  of  observations  are  made  on  separate  flights, 
the  navigation  vectors  Sj  and  Sg  are  measured,  and  S  involves  the 
errors  in  both  vectors. 

For  the  purposes  of  an  error  analysis,  the  complete  solution  to 
Equation  4-1  is  not  required.  The  error  in  P  due  to  an  error  in  a  meas¬ 
urement  of  variable,  may  be  found  by  differentiation.  If  ^  is  the  <  th 
omponent  of  vector  ^  and  /Uj  is  the  j  th  measured  variable  then; 

4-4 

where  tbe  error  in  due  to  the  error  in  /l/J  The  total 

error  in  pj  can  be  found  from  (dPi)j  Assuming  that  the are 
chosen  to  be  directly  measured  quantities,  the  errors  in  different  are 

independent,  Statistically,  if  the  in  Equation  4-4  is  the  probable  error 

in  ,  then  the  probable  error  in  ^  due  to  all  error  sources  is  given 
by: 


If  an  over-all  measure  of  target  location  accuracy  is  desired,  the 
root  mean  square  error  in  the  magnitude  of  the  distance  between  the  actual 
target  and  its  measured  location  may  be  calculated.  Consideration  of  vector 
error  shows  this  to  be  simply: 


Combining  Equations  4-4,  4-5  and  4-6: 

Therefore,  the  required  quantity  is  the  ^  matrix.  The  vector  it 
given  by  Equation  4-1,  which  may  be  differentiated  with  respect  to  ,  cox 
ponent  by  component.  This  process  yields: 


In  some  cases,  the  horizontal  or  vertical  error  in  target  location  may  be  of 
interest.  These  separate  errors  may  be  calculated  by  changing  Equation  4-6, 
All  other  parts  of  the  analysis  are  applicable. 
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4-8 


b/LTf 


\  . 

Since  the  components  of  S)  are  directly  measured,  unless 

-CS,)f  >  and  s  /  .  The  matrix  must  be  evaluated  for 

each  sensor  combination,  depending  on  the  particular  information  about 


and  fSg  which  is  measured  by  those  sensors. 


For  the  parametric  study  described  later,  three  sensor  combinations 
are  used.  These  are:  Panoramic  camera  stereo,  side-looking  radar  stereo, 
and  a  panoramic  camera  combined  with  side-looking  radar.  The  derivations 
of  ''  values  for  these  cases  are  given  in  Appendices  A,  B,  and  C. 

Other  possible  sensor  combinations  would  include  Frame  cameraand  IR. 

The  long  focal  length  cameras  considered  restrict  frame  cameras  to  narrow 
field  angles.  Over  these  angles,  the  geometry  in  frame  and  panoramic  cameras 
are  nearly  identical,  since  the  curved  focal  surface  only  affects  larger  angles. 
Therefore,  the  panoramic  camera  results  at  small  angles  may  be  assumed  to 
be  similar  to  what  would  be  obtained  for  frame  cameras.  IR  was  not  con¬ 
sidered  since  the  present  state-of-the-art  does  not  provide  for  high  altitude 
capability.  However,  including  IR  results  requires  a  straightforward 
extension  of  the  calculations. 


•  1 

.V\i 


5.  PARAMETRIC  STUDY  OF  TARGET  LOCATION  ERRORS 


5.  1  General  Description 

The  target  location  study  was  not  to  be  tied  to  any  particular 
state-of-the-art  of  any  equipment  involved,  so  the  measurement  errors  were 
varied  from  presently  expected  values  down  to  zero,  to  allow  for  future 
improvements.  Tho  altitude  ranges  considered  are  the  regime  of  boost- 
glide  type  vehicles.  The  altitude  range  limits  were  50,000  and  1,500,000 
feet. (Thus  including  high  altitude  cruise  type  vehicles,  as  well  as  satellites.) 

This  parametric  treatment  of  the  problem  has  several  purposes. 
Primarily,  it  is  desired  that  the  trade-offs  between  various  combinations  of 
sensors  could  be  illustrated  as  a  function  of  measurement  erro.  s  and  param¬ 
eters  such  as  sensor  altitude.  These  comparisons  and  trade-offs  are 
treated  after  the  description  of  probable  error  calculations. 

It  is  also  possible  to  use  the  output  of  this  study  to  predict  the 
probable  error  in  targot  location  when  any  particular  sensor  combination 
is  used.  The  error  introduced  by  various  measurement  errors  is  provided 
as  a  function  of  geometrical  factors  (altitude,  etc.  )  for  each  sensor  combination. 

From  the  target  location  error  results,  it  is  also  possible  to  draw 
conclusions  with  respect  to  image  restitution  for  purposes  of  image  regis¬ 
tration. 

The  parametric  study  described  in  the  remainder  of  this  section 
was  carried  out  using  the  Laboratory's  IBM  704  computer.  The  ‘target 
location  errors  are  computed  by  making  use  of  the  location  error  equations 
developed  in  the  previous  section  and  its  appendices.  The  results  are  in  the 
form  of  graphs  for  ease  of  examination.  Two  sets  of  graphs  are  presented; 
one  which  shows  the  contribution  to  target  location  error  due  to  various 
measurement  errors,  as  a  function  of  location  geometry,  and  a  set  which 
shows  trade-offs  between  systems  as  a  function  of  sensor  measurement 
errors. 
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5.  2  Error  Contributions  due  to  Various  Sensor  Errors, 


Three  basic  sensor  combinations  were  chosen  for  the  study.  These 
are:  Stereo  Panoramic  camera,  Panoramic  camera  with  side-looking  radar, 
and  Stereo  radar.  The  frame  camera,  an  obvious  possible  sensor  choice, 
has  not  been  included,  since  the  geometry  of  frame  and  panoramic  systems 
is  sufficiently  similar  at  small  field  angles  that  the  error  analysis  is  not 
affected  by  the  choice.  The  panoramic  camera  was  chosen  rather  than  the 
frame  camera  in  order  to  Include  larger  field  angles  as  well  as  the  small 
angles  of  a  frame  system.  Both  front-  and  side-lap  is  considered  for  stereo 
measurements  where  applicable.  (High  resolution  radar,  at  present,  cannot 
be  used  for  stereo  by  measuring  with  forward  and  rearward  squinted  antennas 
on  a  single  flight  line.  ) 

Throughout  the  derivations  of  target  location  equations  in  Section  4 
camera  errors  were  assumed  to  be  the  errors  in  film  measurement  and  focal 
length.  It  is  clear  that  these  are  not  the  only  sources  of  error.  However, 
other  errors  can  be  included,  since  other  errors  can  only  affect  the  direction 
of  the  camera  vector.  Just  as  film  and  focal  length  errors  do.  For  example, 
atmospheric  hydrostatic  deflection  effectively  distorts  the  image,  and  is  there¬ 
fore  equivalent  to  a  film  measurement  error.  Thus  the  errors  in  film  and 
focal  length  measurement  are  sufficient  to  take  into  account  all  possible 
camera  errors. 

In  order  to  reduce  the  number  of  separate  error  sources  still 
further,  so  that  the  total  output  of  information  is  not  overwhelming,  other 
error  sources  were  grouped  together.  All  navigation  errors  (error  in  S,  S, 
and  Sx  of  Section  4}  were  grouped  and  called  simply  navigation  errors.  Like¬ 
wise,  all  roll,  pitch  and  yaw  errors  which  contribute  to  location  error  are 
included  as  stabilization  errors.  Film  and  focal  length  errors  include  those 
errors  due  to  all  cameras  in  the  system.  In  the  case  of  radar  stereo,  the 
antenna  beam  spread  in  the  azimuthal  direction  adds  directly  to  yaw  errors 
in  target  location,  so  it  is  included  as  a  stabilization  error.  The  total 
number  of  groups  of  errors  is  therefore  reduced  to  5,  namely: 

1}  Radar  range  error 
2}  Film  measurement  error 
3)  Focal  length  error 
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4)  Navigation  error 

5)  Stabilization  error 


These  five  errors  can  include  all  possible  error  sources  in  a  system.  Table 
5.  1  lists  several  examples. 

The  total  target  location  error  can  be  found  from  the  error  con¬ 
tributions  by  a  simple  procedure.  Each  curve  plotted  in  Figures  5.  3  to 
5.  20  shows  the  error  contribution  due  to  one  of  the  5  errors  listed  above. 

The  error  contribution  is  computed  for  unit  error  in  each  case,  but  for  con¬ 
venience  in  plotting,  the  values  are  then  multiplied  by  constants,  These 
constant  multipliers  are  equivalent  to  assuming  errors  as: 

1)  Range  measurement  =  ^  10  meters 

2)  Film  measurement  =  -f  50  microns 

3)  Focal  length  measurement  =  +  100  microns  5-2 

4)  Navigation  component  =  _+  10  meters 

5)  Stabilization  angles  =  ^  0.  1  milliradians 

If  one  wished  to  provide  different  estimates  of  these  errors,  he 
would  compute  total  location  error  by  reading  the  value  of  each  error  con¬ 
tribution  from  the  appropriate  graph,  dividing  by  the  constant  listed  in  5-2 
to  get  the  contribution  from  unit  error,  and  multiplying  by  his  estimate  of 
error  for  that  measurement.  This  yields  the  location  error  due  to  each 
error  source.  The  total  error  is  obtained  by  taking  the  square  root  of  the 
sum  of  the  squares  of  individual  errors. 

Defining: 

dF,  dN,  dSTAR  =  errors  1  to  5  listed  in  5-  1 ,  measured  in 

units  given  with  values  in  5-2. 

^,1  Cg  ,  ,  Cy  =  error  contributions  due  to  errors  1  to  5, 

the  values  being  read  from  the  graphs. 


Then,  the  target  location  error  is  given  (assuming  normal  distribution  of  errors) 
by: 

./oo. 


Error  (in  meters) 


rrc,  d/?]  J  rc,  dx'^  TCgd  rc^  a  a/|  *  dsTA^ 

IL /o  J  [d-o\  [/oo\i/o\[  a/  . 


5-3 
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Table  5.  1  TYPICAL  ERROR  SOURCES 


1.  Radar  Range  Error 

1.  Timing  error 

2.  Finite  pulse  length  effect 

3.  Processing  errors 

Camera  Errors 

2.  Film  Measurement  -  errors  which  vary  over  the  format 

1.  Film  measurement 

2.  Nonuniform  film  shrinkage 

3.  Lens  distortion 

4.  Atmospheric  distortion 

hydrostatic  gradient 

refraction  at  shock  wave,  boundary  layer 
and  intermediate  region 

5.  Window  distortion 

6.  Variations  of  scan  angle  such  as  caused  by  off-center 
mounting  of  scanning  prism  and  other  mechanical 
errors. 

7.  Improper  film  flattening 

3.  Focal  Length  -  errors  in  scale  which  do  not  vary  over  the 
format,  but  may  vary  from  picture  to  picture 

1.  Focal  length  determination 

calibration 
tolerances  in  camera 

expansion  caused  by  temperature  changes 

2.  Uniform  film  shrinkage 

4.  Navigation 

1.  Errors  in  inertial  system 

2.  Errors  in  marking  of  film  or  tape 

3.  Errors  in  control  point  locations 

5.  Stabilization 

1.  Errors  in  stabilizing  mount 

2.  Errors  in  marking 

3.  Radar  synthetic  antenna  beam  width 
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An  example  of  this  process  is  given  in  Figures  5.  1  and  5.  2. 

The  error  contributions  are  plotted  as  functions  of  five  parameters. 
These  are;  the  altitude  of  sensor  1  ,  the  altitude  of  sensor  2  (^z), 

the  field  angle  or  angle  from  the  vertical  to  the  line  between  sensor  1  and 
the  target  (A),  the  camera  focal  length  (F),  and  the  side  offset  cf  sensor  2 
with  respect  to  sensor  station  1  (5^)  ,  The  last  parameter  (Sy)  is  used  only 
as  zero  or  nonzero  to  separate  forward  and  side-lap  location. 

It  is  obviously  useless  to  plot  all  combinations  of  the  five  param¬ 
eter  values.  Instead,  a  single  set  of  values  was  picked  as  a  starting  point, 
and  all  parameters  were  varied  separately,  holding  the  others  at  the  chosen 
values.  The  chosen  set  of  values  used  as  a  basis  is: 

>/;  =100  Kft 

/Vz  =  100  Kft 

A  =  30*  5-4 

=  1  m 

Sy  -  0  OX  fo  ,  plotted  separately 

The  parameters  are  varied  through  the  following  ranges,  as 
suggested  by  the  sponsoring  agency: 

=  |50  -  1500  Kft  for  cameras 

I  50  -  100  Kft  for  radar  5-5 

A  =  0-60  degrees  from  vertical 

r-'  -  0.  15  •  2  meters 

The  maximum  allowable  baseline  for  forward-lap  photography  was 
set  by  a  choice  of  a  9-inch-wide  format,  which  limits  the  target  angle  ahead  or 
behind  the  camera  which  still  appears  on  the  film.  The  side-looking  radar  was 
assumed  to  have  the  antenna  beam  set  at  right  angles  to  the  flight  direction. 

No  restriction  is  set  on  the  field  angle  at  which  the  radar  can  operate,  but 
it  must  be  realized  that  at  very  small  angles  the  resolution  will  not  be  suf¬ 
ficient  to  identify  a  target.  As  discussed  previously,  it  has  been  assumed 
throughout  the  target  location  analysis  that  the  target  is  identifiable  on  each 
sensor  record. 
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Eimt  COITRIBUTlOa 


ALTITUBC.  k  ft 

Figure  6.1  PAN-RAO  SYSTEM  -  VARIATION  WITH  ALTITUDE  OF  SENSOR  ONE 
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Total  Location  Error  Calculation 


Assume:  H)  -  Hi  -  100  Kft 

A  =  30* 

P  =  1  m. 

Sy  =  0 

et^  -  3  m. 

4/X  =  10 

dP  =  100 
(JN  =  30  m. 

d^TAti  =  0.  1  milliradian 

From  graph  Ct  ~  10.051  m. 

Cl  =  i.489 
C)  -  1.877 
=  17.380 
-  5.016 

From  Equation  5-3: 

Error  ( [(/.877)(/00)]^  r(/73a0)(30)]^  [iS^OI6»){0. 

/ 1  /a  \  [  sro  \  \  too  \  \  to  ]*  o.t 


Error  =  53  m 


Figure  5.  2  EXAMPLE  OF  ERROR  CALCULATION 
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EiTor  Contribution  Results 


5.  3 

5.  3.  1  Panoramic  Camera  Stereo  ("Twopan") 

Considering  the  forward-lap  photographs  first,  it  can  be 
seen  that  as  the  altitude  of  the  cameras  is  increased,  the  error  due  to  any¬ 
thing  except  navigation  is  proportional  to  altitude,  while  navigation  error 
is  altitude  independent,  Figure  5.  3.  Varying  field  angle  has  very  slight  effect 
on  error,  with  a  slight  increase  in  location  errors  at  larger  angles,  Figure  5.4. 

Varying  the  focal  length,  Figure  5.  5,  has  several  effects. 

A  shorter  focal  length  means  smaller  scale,  and  correspondingly  allows  a 
longer  base  line(separation  of  the  two  camera  stationB)for  triangulation.  There¬ 
fore,  navigation  and  stabilization  errors  have  less  effect  on  short  focal  length 
systems.  A  fixed  error  in  focal  length,  of  course,  produces  a  greater  percent 
change  at  short  focal  length,  thus  contributing  more  to  location  error.  The 
effects  of  scale  change  and  base  line  change  effectively  balance  for  film  errors, 
causing  film  error  contribution  to  be  almost  independent  of  focal  length. 


ALTITUDE,  k  ft 

Figure  6.3  TW0~PAN  SYSTEM  -  VARIATION  WITH  ALTITUDE 


22 


CM-1418-G-1 


The  case  of  side-lap  stereo  is  similar  to  forward-lap. 
Because  of  slightly  different  geometry,  the  curves  are  not  identical,  but  the 
same  general  character  ii'  shown,  Figures  S.  6  and  5.  7.  In  the  case  of 
variation  of  focal  length,  the  base  line  for  side-lap,  Figure  5.8,  was  kept 
fixed,  so  that  the  curves  for  contribution  due  to  navigation  and  stabilization 
become  almost  independent  of  focal  length,  while  film  error  contribution 
follows  the  same  form  as  focal  length  contribution,  due  to  scale  change. 


ALTITUDE,  li 

Figure  6.6  TWO-PAH  SYSTEM  VARIATION  WITH  ALTITUDE 
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FIgurt  6.7  T«0>FAN  SYSTEM  -  VARIATION  WITH  FIELD  ANGLE 


Figure  5.8  TWO-PAN  SYSTEM  -  VARIATION  WITH  FOCAL  LENGTH 


j 


5 
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5.3.2  Panoramic  Camera  with  Radar  ("Panrad") 


A  system  combining  a  camera  with  radar  makes  use  of  the 
direction  determination  by  the  camera  and  the  range  determination  of  the 
radar.  If  the  photograph  is  taken  from  nearly  the  same  place  as  the  radar 
record,  the  combination  provides  range  from  the  radar  and  angles  from  the 
camera  to  fix  a  single  vector.  If  the  two  sensors  are  used  far  apart,  the 
situation  is  more  complicated.  The  target  location  error  becomes  extremely 
large  if  the  lines  from  the  target  to  each  sensor  meet  at  right  angles,  since 
the  location  is  then  found  as  the  intersection  of  a  sphere  and  a  tangent  line. 
Thus,  for  certain  base  lines,  the  system  fails,  Just  as  photographic  stereo 
fails  if  the  base  line  is  reduced  toward  zero. 

As  in  the  photographic  stereo  system,  then,  the  two  sensor 
stations  can  be  either  essentially  on  the  same  flight  path  (SysOJor  on  separated 
paths  (5/ ^  C)  ,  As  previously  discussed,  the  radar  sensor  altitude  was  held 
at  100  Kft  or  lees.  This  means  that  the  cases  listed  as  Sy*0are  situations 
where  the  camera  and  radar  sensors  are  carried  either  at  the  same  altitude 
(possibly  1  vehicle)  or  else  that  the  two  flights  are  at  different  altitudes,  but 
are  close  enough  to  each  other  horizontally  that  the  target  is  on  the  same  side 
of  both  flight  lines. 

Examining  Figure  5.9,  it  can  be  seen  that  for  5y0  the 
error  contributions  as  a  function  of  camera  altitude  are  much  like  those  for 
two  camera  stereo,  with  the  addition  that  radar  range  error  contribution  is 
essentially  independent  of  camera  altitude.  Figure  5.  10  shows  that  target 
location  error  is  independent  of  radar  sensor  altitude  over  a  factor  of  two  in 
altitude. 
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As  with  2  cameras,  a  camera  combined  with  a  radar  on 
nearly  the  same  flight  path  (S/«0)  yields  target  location  errors  which  are 
nearly  independent  of  field  angle.  It  is  apparent  from  Figure  5.  11,  however, 
that  the  contribution  to  location  error  due  to  focal  length  error  increases  with 
field  angle,  becoming  more  critical  at  larger  angles. 

The  system  requirements  for  camera  error  limits  become 
less  severe  as  the  photographic, scale  is  increased.  This  can  be  seen  in 
Figure  5.  12  showing  error  contribution  as  a  function  of  focal  length  of  the 
camera.  As  focal  length  is  increased,  the  scale  becomes  larger,  and  the 
location  error  due  to  camera  measurements  becomes  less. 
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WITH  EimOH  VALUES  FROM  TABLE  S.2 
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Flflurt  6.11  PAN-RAD  $mEM  -  VARIATION  WITH  FIELD  ANGLE 
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Figure  8.12  PAN-RAD  SYSTEM  -  VARIATION  WITH  FOCAL  LENGTH 
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If  the  camera  and  radar  are  carried  on  separate  flights 
which  pass  opposite  sides  of  the  target  (Sy/  O),  the  poor  location  geometry 
that  may  occur  can  cause  large  errors. 

As  camera  altitude  is  varied,  Figure  5.  13,  the  same  general 
trends  occur  as  with  Sy  •  0  ,  but  the  poorer  geometry  causes  over-all  larger 
error  contributions.  Varying  the  radar  altitude,  keeping  the  target  location 
with  respect  to  the  camera  fixed,  Figure  5.  11,  causes  the  two  position 
vectors  to  intersect  at  90*  when  the  radar  is  at  about  76  Kft.  Thus  the  error 
becomes  extremely  large  at  that  point. 
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WITH  ERIIOR  VALUER  FROM  TABLE  S.2  t 
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Figuru  5.14  PAN-RAO  SYSTEM  -  VARIATION  WITH  ALTITUDE  OF  SENSOR  TWO 
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Similarly,  as  the  field  angle  is  varied.  Figure  5.  15,  the 
critical  geometry  is  reached  at  about  37*.  again  causing  extremely  large 
errors.  This  corresponds  to  the  poor  geometric  situation  which  arises  when 
two  photographs  to  be  used  for  stereo  are  taken  from  nearly  the  same  point. 

Regardless  of  geometry,  error  contributions  due  to  radar  range, 
camera  stabilisation  and  navigation  errors  are  independent  of  scale,  and  camera 
errors  contribute  less  as  the  scale  is  increased  by  Increasing  focal  length, 

Figure  5.  16. 
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Figure  5.16  PAN-RAO  SYSTEM  -  VARIATION  WITH  FOCAL  LENGTH 
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5.3.3  Radar  Stereo  ("Tworad") 


Radar  stereo  relies  upon  a  knowledge  of  two  ranges  to  locate 
a  target  on  a  circle;  the  target  is  then  located  since  it  must  be  within  the 
radar  antenna  beam  as  well  as  on  the  circle.  The  two  radar  sensor 
stations  must  be  on  separate  flights  since  the  side-looking  radar  records  the 
target  image  at  essentially  90*  to  the  flight  line. 

Figure  5.  17  shows  error  contributions  as  a  function  of  sensor 
altitude,  when  the  two  sensor  stations  are  separated  horizontally  (parallel 
flights).  The  only  variation  observed  is  that  due  to  stabilization.  That  error 
becomes  more  critical  at  higher  altitude  only  because  the  target  field  angle 
was  kept  fixed,  so  the  target  was  farther  away  horizontally  when  the  sensor 
was  higher.  Since  the  radar  antenna  beam  is  used  to  help  locate  the  target,  a 
yaw  error  becomes  more  critical  for  greater  horizontal  distance  to  the  target. 

Figure  5.  18  shows  a  similar  effect  for  stabilization  error 
as  the  horizontal  distance  is  increased  by  increasing  field  angle.  Since  side¬ 
looking  radar  will  not  work  at  small  angles  from  vertical,  the  navigation  and 
range  error  contributions  become  large  at  small  angles. 

If  an  attempt  is  made  to  use  two  radars  at  different  altitudes, 
but  nearly  the  same  horizontal  position,  location  errors  increase  as  vertical 
separation  decreases,  as  shown  in  Figure  5.  19, 

Table  5.  Z  NOMINAL  ERROR  VALUES 

1)  Range  measurement  =  ■)■  10  meters 

Z)  Film  measurement  =  +  50  microns 

3)  Focal  length  measurement  =  ^  100  microns 

4)  Navigation  component  =  +  10  meters 

5)  Stabilization  angles  =  ■•■  0. 1  milliradians 
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ALTITUDE,  k  ft 

Figure  5.19  TWO-RAD  SYSTEM  •  VARIATION  WITH  ALTITUDE  OF  SENSOR  TWO 
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5.4 


Comparisons  of  Different  Sensor  Combinations 

In  the  previous  section,  the  results  of  the  parametric  study  were 
presented.  However,  the  form  in  which  the  results  are  presented  does  not 
lend  itself  to  easy  comparisons  of  different  system  capabilities,  Therefore, 
a  few  representative  situations  were  chosen,  and  the  total  target  location 
error  was  plotted  as  a  function  of  various  system  errors,  allowing  trade¬ 
offs  between  sensor  combinations  to  be  shown  directly. 

The  sample  cases  picked  were:  high  and  low  altitude  cameras 
used  for  forward  and  side-lap  stereo,  and  corresponding  cases  for  the  other 
two  sensor  combinations  used.  Again,  radar  is  limited  to  operate  at  100  Kft 
or  less.  Four  cases  are  plotted  for  each  sensor  combination.  The  param¬ 
eter  values  are  given  in  Table  5.  3. 

In  order  to  further  reduce  the  large  number  of  cross-plots, 
system  errors  were  combined  in  two  more  cases.  First,  since  the  only 
stabilisation  error  which  enters  for  radar  stereo  is  yaw  error,  and  since  even 
using  coherent  radar  the  effective  antenna  beam  width  is  the  order  of  a  milli- 
radian  or  more,  the  combined  error  due  to  yaw  and  antenna  beam  spread  was 
held  constant  at  1  milliradian  for  the  study. 

Second,  to  obtain  a  single  number  for  camera  error,  define; 

ax  s  error  contribution  due  to  error  (dXf)  in  film  measurement 

a.pcLf’  s  error  contribution  due  to  error  (dF  )  in  focal  length  meas¬ 

urement 

Since  these  errors  are  independent,  they  should  be  added  as  root  mean  square 
values.  The  result  is: 


/e.e.  ‘  [a/dx^  5-6 

But  the  last  equation  may  be  rewritten: 
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Since  film  measurement  error  is  usually  small  compared  to  focal  length  error, 
the  term  in  brackets  is  not  much  different  from  unity  in  many  cases.  Camera 
error  is  defined  as: 


Each  graph  is  plotted  showing  target  location  error  versus  camera  error  (as 
Just  defined)  for  each  system,  with  radar  error  as  a  parameter.  Separate 
graphs  show  cases  1  to  IV  and  different  navigation  and  stabilization  errors. 

Errors  were  chosen  in  the  ranges: 

Camera  error  •  0  to  lOOOyU 

Radar  error  >  (  1  to  10  meter  range  error 

(  1  milliradian  beamwidth  5-9 

Navigation  error  -  1  to  30  meters 
Stabilization  error  •  .01  to  .  1  milliradian 

Navigation  errors  can  well  exceed  30  meters,  of  course,  but  the  results  show 
that  larger  errors  than  this  will  override  the  effects  of  all  other  errors.  This 
can  also  be  seen  by  referring  to  any  of  the  graphs  in  subsection  5,  3  where  a 
navigation  error  of  10  meters  is  seen  to  contribute  errors  which  are  larger 
than  or  at  least  comparable  to  any  other  error  source. 

5.  5  Results  of  System  Comparisons 

The  results  of  the  system  comparison  study  are  presented  graph¬ 
ically  at  the  end  of  this  section.  The  desirable  system  to  be  used  for  target 
location  is  that  system  which  produces  the  smallest  target  location  error, 
that  is,  the  system  which  is  represented  by  the  lowest  curve  on  a  figure. 

Frequently  two  curves  cross,  in  which  case  the  better  system  is  the  one  for  which 
the  location  error  curve  is  lower  at  the  value  of  camera  error  which  is  attainable. 

Using  Figure  5.  20  as  an  example,  it  can  be  seen  that  stereo  panoramic 
cameras  (2  Pan)  are  the  berft  location  system  as  long  as  camera  error  is  less  than 
720  microns.  If  camera  error  is  greater  than  that  value  (probably  due  to  focal 
length  errors),  radar  stereo  (2  Rad)  would  be  better,  if  range  error  can  be  held 
to  1  meter.  The  camera- radar  (Panrad)  combination  is  preferred  over  2  Rad 
only  when  camera  error  is  less  than  80  microns  with  radar  range  error  of  1  meter. 
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Figures  5,  21  to  5.  23  show  low  altitude  (100  Kft),  separated, 
parallel  flights  over  a  target  (Case  I,  Section  5.4).  Comparison  of  5.  20  and 
5.  21  show  that  stabilization  error  does  not  appreciably  affect  results  at 
these  attitudes  if  held  under  0.  1  milliradian.  which  should  be  possible. 

Figures  5.  21  to  5.  23  illustrate  the  effect  of  increasing  navigation  error. 

It  is  clear  that  Panrad  is  not  to  be  used  for  this  situation,  and  that 
photography  is  better  than  radar  so  long  as  camera  error  is  less  than  about 
700  microns.  If  navigation  errors  become  larger  than  about  10  meters,  radar 
and  photographic  stereo  yield  comparable  results. 

It  should  be  noted  that  navigation  and  stabilization  errors  need  not 
be  as  large  as  would  be  the  case  for  self-contained  systems.  The  errors  may 
be  reduced  by  using  control  extension  techniques,  and/or  ground  based  tracking 
systems.  These  reduced  navigation  errors  can  be  used  with  the  graphs  pre¬ 
sented  here,  to  evaluate  any  target  location  system. 

Figures  5,  24  to  5.  27  show  the  location  errors  when  one  flight 
line  is  used  (or  one  flight  above  another  for  radar).  It  can  be  seen  in  general 
that  this  geometry  is  ideal  for  a  panoramic -radar  combination  and  that  this 
system  should  outperform  even  a  2  photo  determination,  provided  that  equi¬ 
valent  navigational  and  orientation  data  are  available.  The  photographic 
stereo  will  only  produce  better  results  if  camera  error  is  very  small  and 
navigation  data  is  excellent,  If  navigation  errors  exceed  10  meters,  radar 
becomes  preferable  to  photographic  stereo,  since  radar  target  location  is 
less  affected  by  navigation  errors. 

The  cases  in  which  the  cameras  are  carried  at  high  altitude 
(1,000  Kft)  are  shown  in  Figures  5.28  to  5.37. 

The  discussion  of  the  low  altitude  cases  can  be  applied  as  well 
to  high  altitudes,  except  that  given  navigation  and  stabilization  errors  intro¬ 
duce  larger  target  location  errors  than  at  lower  altitudes. 
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Fiquro  5.20  TARQET  LOCATION  ERROR 
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Figure  5.21  TARQET  LOCATION  ERROR 
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Figure  6.30  TARGET  LOCATION  ERROR 
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Figure  5.31  TARQLT  LOCATION  ERROR 
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Figure  6.32  TARGET  LOCATION  ERROR 


Figure  5.33  TARGET  LOCATION  ERROR 
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6.  CONCLUSIONS 


Although  the  study  is  primarily  concerned  with  target  location 
and  charting,  irttportant  conclusio'ns  about  image  registration  must 
be  included: 

a)  For  t7ne  altitude  range  of  boost-glide  vehicles  considered  here, 
relief  displacement  errors  in  single  images  prevent  registration  of  all  except 
large  targets.  Even  where  stereo  methods  are  used  point  location  errors  are 
still  too  large  for  images  to  be  registered  by  only  rectification  of  the  images. 
Improvement  can  be  achieved  only  by  identifying  several  points  on  each  stereo 
model  and  forcing  a  local  matching. 

Conclusions  concerning  target  location  and  charting: 

b)  Navigation  errors  in  excess  of  about  30  meters  cause  all  systems 
to  be  limited  by  that  error.  If  inherent  sensor  capabilities  are  not  to  be 
overwhelmed  by  navigation  errors,  present  systems  require  the  aid  of  ground 
control  data  or  control  extension  techniques. 

c)  As  in  the  case  for  photographic  systems,  the  best  geometry  for 
radar  stereo  systems  is  that  which  causes  the  lines  from  the  two  sensors 
to  the  target  to  meet  at  right  angles.  For  a  camera-radar  combination,  the 
best  geometry  occurs  with  the  two  sensors  carried  in  one  vehicle.  If  both 
sensors  are  carried  in  one  vehicle,  such  a  combination  can  outperform  other 
target  location  systems. 

d)  Stabilisation  errors  do  not  appreciably  affect  low  altitude  systems, 
if  present  stabilization  errors  are  assumed.  At  higher  altitudes,  stabilization 
errors  can  dominate  in  photographic  systems. 


Rectification  Includes  all  corrections  which  can  be  made  by  calibration  or 
calculation.  For  example,  scale  variation  over  the  field,  lens  distortion  and 
some  static  atmospheric  effects  are  corrected  during  this  process. 

Stereo  is  used  in  the  broad  sense  defined  at  the  start  nf  Part  I  of  this  report. 


51 


CM-1418-G-1 


e)  When  a  large  baae  line  (separation  of  sensor  stations)  is  avail¬ 
able,  as  when  separate  flights  or'  convergent  photography  is  used,  photographic 
stereo  provides  the  most  accurate  target  location. 

f)  Since  radar  is  essentially  independent  of  stabilization,  and  tri- 
lateration  with  radar  stereo  is  only  slightly  dependent  on  relative  positions 
of  the  two  sensors,  a  trade-off  point  exists  between  radar  and  photographic 
stereo  systems.  The  radar  system  outperforms  the  photographic  system 
if  navigation  and  stabilization. errors  are  large. 
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7.  RECOMMENDATIONS 


The  present  study  provides  a  comparison  of  the  capabilities  of 
sensor  combinations  for  given  values  of  error  due  to  various  sources.  This 
study  evaluates  target  location  error  in  terms  of  basic  target  location  equations. 

Image  registration  by  rectification  only  has  been  shown  to  be 
impractical,  even  if  stereo  models  are  prepared  first.  Two  other  possibilities 
should  be  studied.  First,  the  improvement  possible  by  identifying  points  on 
each  model  and  forcing  these  corresponding  points  to  register  should  be 
investigated.  Such  local  sc  ale -changing  might  register  the  images  quite  well 
over  the  area  surrounding  the  identifiable  points.  Secondly,  other  possible 
uses  of  less  precise  registration  could  be  studied.  For  example,  radar  could 
be  used  for  target  detection,  and  a  coarse  image  registration  would  allow  a 
small  area  of  a  photograph  to  be  isolated,  with  the  target  known  to  be  within 
this  area. 

Research  to  improve  navigation  systems  should  be  continued,  since 
present  navigation  errors  limit  sensor  capabilities. 

Further  study  should  be  carried  out  to  determine  the  precise  reduc¬ 
tion  of  target  location  error  when  redundant  information  is  Included  as  well  as 
when  control  extension  techniques  are  used.  Such  study  could  follow  the  outline 
of  the  least  squares  analysis  presented  in  this  report,  but  should  include  control 
extension. 


See  first  footnote,  p,  51 
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PART  11 


8.  ATMOSPHERIC  EFFECTS 


The  probable  contributions  to  target  location  errors  due  to  atmospheric 
effects  can  be  treated  analytically  and  evaluated  numerically  if  a  few  con¬ 
stants  are  obtained  from  experimental  data.  Steady  refraction  effects  need 
not  contribute  to  location  errors,  since  they  can  be  evaluated  and  taken  into 
account  in  correcting  photographic  measurements  for  systematic  errors 
(rectification).  Fluctuations  in  atmospheric  refraction,  and  random  errors 
in  evaluating  steady  effects  are  accounted  for  in  the  target  location  error 
analysis  in  Part  I.  In  the  following,  all  refraction  effects  are  treated  ana¬ 
lytically,  and  numerical  estimates  are  made  on  the  basis  of  available  data. 

Refraction  effects  in  airborne  photographic  systems  result  from  the 
following  density  gradients  in  the  earth's  atmosphere  and  in  the  flow  field 
associated  with  the  airborne  vehicle: 

1.  The  hydrostatic  density  gradient  in  the  ambient  atmosphere 
due  to  the  earth's  gravitational  field; 

2.  Heating  effect  gradients  in  the  ambient  atmosphere  or  in  the 
flow  field; 

3.  The  density  Jump  across  the  shock  wave  accompanying  a  super¬ 
sonic  aircraft; 

4.  The  density  gradient  associated  with  the  steady,  continuous 
pressure  distribution  in  the  flow  field  related  to  the  vehicle 
shape,  attitude,  and  flight  conditions; 

5.  Random  density  fluctuations  associated  with  turbulence  in 
the  boundary  layer  or  in  the  ambient  atmosphere. 

The  resulting  degradation  of  the  photographic  image  can  be  described 
in  terms  of  the  following  phenomena: 

1.  Displacement  of  the  image  by  steady,  uniform  deflection  of 
all  the  image  forming  rays; 

2.  Warping  of  the  image  by  steady  but  non-uniform  deflection  of 
the  image  forming  rays; 


55 


CM-14i8-G-l 


3.  Shimmy  or  random  motion  of  the  image  points  by  unsteady  but 
uniform  deflection  of  the  rays; 

4.  Distention  or  image  point  spread  by  unsteady  and  non-uniform 
deflection  of  the  ray;; 

5.  Twinkling  or  brightness  flnctuations  due  to  random  deflection 
at  a  considerable  distance  from  the  camera; 

6.  Dispersion  or  separation  of  colors,  either  steady  or  random, 
due  to  the  spectral  dependence  of  the  specific  refractivity. 

The  various  refraction  effects  will  be  treated  in  the  order  listed,  starting 
with  the  hydrostatic  density  gradient. 

8.  1  Hydrostatic  Gradient 

The  relation  between  the  variation  of  the  index  of  refraction  and 
the  resulting  deviation  is  given  by  the  curvature  equation,  Reference  1, 

where  i*  is  the  direction  vector  tangent  at  any  point  to  the  image  forming 
light  ray  passing  through  that  point  and  is  the  curvature  of  the  ray. 

For  most  purposes  the  index  of  refraction  n  can  be  expressed  in  terms  of 
a  constant  specific  refractivity/(|R«feronces  2,  3i  4,5, 

n  I  Kp 

where  8-2 

.  23  cm^/gm  b  .  12  ft^/slug 

For  aerial  photography  of  the  ground  from  altitudes  below  300,000 
feet  the  direction  of  the  vertical  varies  less  than  10*  along  the  ray  from  any 
visible  point  on  the  ground  to  the  aircraft  so  that,  for  practical  purposes,  the 
earth  can  be  treated  as  flat,  the  hydrostatic  density  gradient  as  unidirectional, 
and  Snell's  law  can  be  used  to  describe  the  angular  displacement  6  at  the 
airplane  of  the  corresponding  level  displacement  d,  along  the  ground  by  the 
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equation 


eh 

'ta,r\  aC(, 


M/ 


8-3 


where  px  iB  the  ambient  density  at  the  aircraft  altitude,  the  summation  is  over 
the  layers  of  constant  lapse  rate  Lih)  from  the  ground  to  the  aircraft,  A/  are 
the  coefficients 


A. .  -JUi 

and  Ij  and  are  the  integrals 


8-4 


if  L‘  O 

8-5 

if  L  ^0 


and 


8-6 


The  reference  values  and  coefficients  for  the  ARDC  Model  Atmosphere,  1959 
are  given  in  Table  8.  1  and  the  generalized  displacement  expressed  in  Equation 
8-3  is  plotted  in  Figure  8.  1. 
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h  (hft) 
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TABLE  8.  1 


THE  ARDC  MODEL  ATMOSPHERE,  1959 
(Reference  6) 


Hi' 

L*/0‘ 

Ti 

ft 

9 

C| 

A  «■ 

(Kft) 

(•R/Kft) 

CR) 

(slug/ft^)  (Ib/slug) 

( 'R/Kft) 

0 

3.  59 

519 

2.  377  X  10"^ 

32.  17 

15.  13 

.  237 

17350 

36 

0 

390 

7.  055  X  10"'* 

32.06 

18.  62 

0 

2520 

82 

-1.63 

390 

7.  761  X  10'® 

31.92 

20.  16 

-.  081 

1080 

155 

0 

509 

2.  883  X  10*^ 

31.70 

18.43 

0 

3310 

175 

2.  85 

509 

1. 396  X  10*^ 

31.  64 

15.  52 

.  184 

46700 

249 

0 

298 

8.  131  X  10-® 

31.42 

18.  23 

0 

1960 

300 

298 

4.  256  X  10-9 

31.  27 

Deflection  of  the  light  rays  by  the  hydrostatic  density  gradient  is 
steady  and  results  in  a  radial  warping  about  the  nadir  point.  Since  the  displace¬ 
ment  is  positive,  targets  will  be  closer  to  the  nadir  point  than  indicated  by 
the  straight  line  tangent  to  the  light  ray  at  the  camera  or  the  image  is  expanded 


radially  from  the  image  of  the  nadir 
point. 
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8.  2  Temperature  Effects 

The  boundary  layer  around  an  aircraft  and  thermal  convection  cur¬ 
rents  in  the  atmosphere  are  regions  of  nonuniform  temperature  but  nearly 
uniform  pressure  where  the  density  varies  inversely  with  the  temperature. 

The  relatively  steady  temperature  gradients  in  the  atmosphere  are  small  and 
have  little  effect  on  light  rays  while  very  high  temperatures  are  reached  at 
the  skin  of  hypersonic  vehicles  and  the  resulting  variation  in  refractivity  may 
significantly  deflect  rays  reaching  the  boundary  layer  with  a  high  angle  of 
incidence. 

At  the  high  supersonic  velocities 
for  which  this  deflection  is  most  appre¬ 
ciable  the  temperatures  at  the  wall 
becomes  so  high  and  the  density  so  low 
that  the  refractivity  n-/  at  the  wall 
vanishes  and  the  variation  of  refractivity 
across  the  boundary  layer  is  nearly  to 
its  value  outside  the  boundary  layer. 

Then 


€  ■  d  •  •ion  (if  8-7 

The  effect  of  this  deflection  on  the  photographic  image  is  a  com¬ 
bination  of  displacement  and  warping  depending  on  the  range  of  the  angle  of 
incidence  fSi  over  the  field  of  view.  The  deflection  and  displacement  are 
positive  indicating  a  stretching  of  the  image.  For  turbulent  boundary  layers, 
the  random  effects  of  shimmy  and  distention  can  be  treated  as  superimposed 
independent  effects.  These  will  be  treated  later. 


8.3  Shock  Wave 


For  the  case  of  a  shock  wave  the  upstream  and  downstream  angles 
of  incidence,  (Sf  and  (Sg,  respectively,  are  related  by  the  expression 


»ii1  g  JIl  ~ 
stnf&f  f)% 


8-8 


60 


CM-1418-G-1 


I 


where  the  deflection  modulus  D  is  defined  by  the  expression 

O-V,  8-9 

and  where  second  order  terms  in  the  refractivity  are  neglected. 

But  sin  can  be  expressed  as  a  function  of  and  the  deflection 

by  the  series  expansion 

stnfSt,  a  smfS,  a  ccsfO/  -  6  — — —  -  &  —  -f  •  •  ■ 

Then,  neglecting  second  and  higher  order  terms  in  the  small  angle  £  we  have 


s/n  fii 


(  *  £  co& 


from  which 


£  •  -D  iAn  fl, 


8-10 


The  angle  of  incidence  /d/  depends  on  the  viewing  angle  and  the 
orientation  of  the  shock  wave.  The  normal  to  the  shock  and  the  upstream  and 
downstream  light  rays  are  coplanar  and  it  is  in  this  plane  that  the  angles  of 
incidence  and  the  deflection  are  measured.  The  shock  wave  angle,  however, 
is  measured  in  the  plane  of  the  normal  to  the  shock  and  the  aircraft  velocity 
vector.  In  the  simplest  case,  these  two  planes  coincide  and  the  angle  of 
incidence  (9^  is  related  to  the  shock  wave  angle  6w  with  respect  to  the 

horieontal  and  the  viewing  angle  of 
with  respect  to  the  vertical  by  the 
expression  •  l&t\/  ~OLf 
where  the  view  angle  aC  is  positive 
for  viewing  forward,  negative  for  rear¬ 
ward,  while  the  angle  of  incidence 
and  deflection  modulus  D  are  positive 
and,  consequently,  the  angular  deflection 
£  is  always  negative,  i.  e. ,  the  upstream  angle  of  incidence  is  greater  than  the 
downstream  and  the  image  is  compressed.  The  geometry  for  shockwaves  of 
general  orientation  is  more  involved  but  the  calculations  are  straightforward 
and  the  details  are  omitted  here. 
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At  low  supersonic  Mach  numbers,  the  shock  wave  does  not  lie  back 
along  the  body  but  is  nearly  perpendicular  to  the  line  of  flight.  In  this  case, 
some  rays  may  cross  the  shock  in  the  manner  already  described,  some  of  the 
rays  from  the  rear  may  reach  the  camera  without  crossing  the  shock,  and 
some  rays  may  be  reflected  to  the  camera  from  the  downstream  side  of  the 
shock.  A  further  conxf^cation  is  the  presence  of  the  trailing  shock  wave  caused 
by  the  convergence  of  the  air  behind  the  vehicle.  Expressions  are  available, 
Reference  7,  for  finding  the  strength,  shape,  and  distance  between  these  shock 
waves  but  the  results  are  mainly  qualitative  and  more  useful  in  determining 
the  optimum  camera  placement  and  flight  speed  rather  than  in  predicting  the 
position  and  extent  of  the  distortion  bands  on  the  photograph. 

8.4  Flow  Field 

Between  the  body  of  a  supersonic  aircraft  and  its  shock  wave,  there 
is  a  region  in  which  the  pressure,  density,  and  fluid  velocity  vary  smoothly 
and  continuously.  The  distribution  of  these  flow  quantities  throughout  the 
field  depends  on  the  shape  and  velocity  of  the  aircraft  and  on  the  ambient 
atmospheric  conditions.  A  discussion  of  the  techniques  for  calculating  the 
density  distribution  in  such  a  field  is  too  involved  to  be  included  hero  except  to 
mention,  as  an  example,  the  method  of  characteristics.  This  method  is  one 
of  the  moat  accurate  and  useful  approaches  to  finding  the  field  quantities 
for -wiper  sonic,  two-dimensional  or  axially  symmetric  flow  and  three- 
dimensional  effects  can  be  included  as  superimposed  cross  flow  effects.  The 
solution  consists  of  two  families  of  intersecting  characteristic  curves  and  the 
values  of  the  flow  quantities  at  the  intersections.  If  these  curves  are  plotted 
and  the  value  of  the  density  given  at  each  intersection,  curves  of  constant 
density  can  be  faired  across  these  curves  and  the  density  gradient  estimated 
on  the  basis  of  the  direction  of  and  distance  between  these  isochors.  Then 
the  deflection  can  be  computed  as  the  product  of  the  specific  refractivity  K  , 
the  average  perpendicular  component  of  the  density  gradient,  and  the  length 
of  the  path.  Again,  this  deflection  occurs  so  close  to  the  camera  that  the 
angular  displacement  can  be  set  equal  to  the  angular  deflection. 


62 


CM-1418-G-1 


8.  5  Random  Fluctuatione  (Turbulence) 


To  treat  random  density  fluctuations  statistically,  we  will  assume 
that  turbulent  regions  in  the  atmosphere  and  the  turbulent  boundary  layer  around 
an  aircraft  are  homogeneous  and  isotropic,  i.  e.  ,  that  the  statistical  quantities 
are  uniform  throughout  and  independent  of  direction.  These  conditions,  of 
course,  are  not  realized  since  the  size,  strength,  and  orientation  of  the  eddies 
depends  partly  on  the  flow  velocity  gradient  which  is  certainly  neither  uniform 
nor  isotropic,  but,  to  a  first  approximation,  the  assumption  is  not  too  unreason¬ 
able. 

A  light  ray  in  a  region  of  turbulent  density  fluctuations  can  be 
pictured  as  a  sinuous  curve  wriggling  around  in  a  very  small  turbular  region 
which  surrounds  the  mean  path  but  which  has  no  sharply  defined  boundary  since 
the  displacement  of  the  ray  from  its  mean  position  is  described  by  a  probability 
distribution.  The  instantaneous  curvature  relation  can  be  expressed  in  the 
form 

where  the  random  curvature  es  well  as  the  deflection 

S'  *  and  fluctuation  n  *  22^  are  random  quantities  while  the 

’  n 

bars  indicate  statistical  means. 

The  random  distribution  of  density  Inhomogeneitles  is  best  discribed 
in  terms  of  a  correlation  function 

8-12 

7* 

where  ^  is  the  distance  between  the  two  points  x'  and  X;  ,  Without 
Justification  except  that  it  has  the  value  1  and  slope  y/^(o)  =  0  at  $  =0, 

decreases  monotonically  to  ^(oc)  =  0,  and  seems  intuitively  to  describe  the 
expected  distribution,  we  shall  use  the  exponential  function 

- a  * 

8-13 

where  the  reference  value  is  /T*  times  the  inflection  value  and  is  char¬ 
acteristic  of  the  turbulent  structure.  With  this  correlation,  the  mean  square 
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deflection  of  a  ray  passing  through  a  homogeneous,  isotropic  turbulent  layer 
of  thickness  s  ,  turbulent  scale  ,  and  mean  square  density  fluctuation 


is  given  by  the  expression 


6^(a)  =  4/C 


m-') 


8-14 


8-15 


The  corresponding  displacement  Q  of  an  image  point  is 
r,'  •  at  6 

where  s/  is  the  image  distance  which,  by  the  lens  law,  is  related  to  the  object 
distance  it  and  the  focal  length  4  by  the  equation 

s,f 

*• 

and  6  is  the  angular  displacement  at  the  camera, 

6n 

where  is  the  distance  from  the  camera  to  the  turbulent  layer.  Then 

For  the  turbulent  boundary  layer,  the  ratio  in  parentheses  is  very  nearly  unity 
and  the  angular  displacement  is  equal  to  the  deflection  while,  for  a  turbulent 
layer  in  the  atmosphere,  the  qviantity  in  parentheses  is  nearly  equal  to  the 
ratio  of  the  altitude  of  the  turbulent  layer  hf  to  the  altitude  of  the  camera  h  , 


for  a  turbulent  boundary  layer  8-16 

for  an  atmospheric  boundary  » 

layer 


Because  the  structure  of  turbulence  in  the  atmosphere  is  not 
homogeneous,  it  is  necessary,  at  least  for  high  altitude  observation  aircraft, 
to  consider  the  atmosphere  as  consisting  of  several  turbulent  layers.  The 
mean  square  angular  displacement  is  then  related  to  the  mean  square  deflec¬ 
tions  by  the  equation 

hi  • 
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since  the  mean  cross  terms 


2  h,  hi  vanish.  Then  the  rms  angular 


displacement  is 


s:  h/ 


and  the  rms  image  point  displacement, 


whereV^i  is  the  rms  deflection  due  to  layer  i 
midpoint  of  the  layer. 


8-18 

and  h  is  the  altitude  of  the 


These  expressions  describe  the  deflection  and  the  displacement  in 
the  image  plane  of  a  single  ray.  But,  unless  the  aperture  of  our  camera  lens 
is  small  compared  to  the  characteristic  size  of  the  density  inhomogeneit.ies , 
we  are  concerned  with  refraction  of  a  beam  of  essentially  parallel  rays  rather 
than  a  single  isolated  ray.  Although  the  random  deflection  of  a  single  ray  is 
the  same  whether  there  are  other  rays  present  or  not,  the  instantaneous 
deflection  of  all  the  rays  is  partially  correlated. 

The  correlation  function  describing  the  correlated 

deflection  of  rays  passing  through  a  circular  aperture  of  radius  Ct  is  given 
by  the  expression 


Kt(k*{)\  (k*3)! 


8-19 


In  the  limit  as  the  aperture  radius  approaches  zero  goes  to  unity  and 

the  instantaneous  average  deflection  covariance  is  equal  to  the  mean  square 
deflection  so  that  the  image  of  a  point  source  becomes,  for  an  ideal  system, 
a  point  moving  over  tlie  circular  region  while  for  large  apertures  the  deflec¬ 
tions  of  widely  spaced  rays  are  uncorrelated  and  the  image  of  a  point  tends  to 
become  a  distended  blur.  The  correlation  function  then  expresses  that 

part  of  the  mean  square  deflection  responsible  for  shimmy,  i,  e.  ,  motion  of 
the  image  point,  so  that  is  called  the  shimmy  factor  and  plotted  in 

Figure  8.  2, 
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FIGURE  8.2  SHIMMY  FACTOR  R\ 


Besides  shimmy  and  distention,  there  is  the  effect  of  twinkle  due 
to  refraction  by  density  inhomogeneities  so  far  from  the  camera  that  the  linear 
displacement  of  the  rays  at  the  camera  is  more  significant  than  the  deflection. 
The  relative  magnitudes  of  these  effects  depend  on  the  aperture  as  well  as  the 
relative  distances  from  the  camera  to  the  turbulence  and  the  object  photo¬ 
graphed  and  can  be  described  in  terms  of  the  quantity, 


,  a* 


8-20 


where  oL  is  the  linear  displacement  of  the  ray  at  the  camera. 

Considerable  twinkle  is  indicated  when  is  of  the  same  order 

of  magnitude  or  greater  than  the  radius  cl  ,  As  an  example,  consider  an  rms 
deflection  of  .02  milliradian  due  to  turbulence  at  20,000  feet,  a  focal  length 
of  1  foot,  and  an  aircraft  altitude  of  40,000  feet.  Then  the  resulting  rms  image 
displacement  is 


which  is  generally  negligible  while  the  displacement  of  the  ray  at  the  objective 

W*  •S(h-ht)  »  0.4  H. 


which,  for  most  camera  apertures,  indicates  considerable  twinkle. 


8. 6  Conclusions 

The  relative  Importance  of  the  various  factors  contributing  to  the 
degradation  of  photographic  images  due  to  refraction  effects  in  the  atmosphere 
and  flow  field  is  indicated  by  Figure  8.  3. 

In  these  curves,  the  most  questionable  data  concerns  the  structure 
and  location  of  turbulence  in  the  atmosphere  and  the  eddy  strengths  and  dis¬ 
tribution  in  turbulent  boundary  layer.  The  values  indicated,  1.  e.  •  /‘^  , 

•  t  P.  *  •"‘i  V,.'  •  were  chosen 

because  they  are  representative  and  are  easily  extrapolated  to  suit  other  condi¬ 
tions.  Since  deflection  by  the  shock  wave  is  proportional  to  and 

tan  fS  *  ,  values  were  chosen  for  P^/P,  and  ^  so  that  these  values  can 

be  conveniently  extrapolated.  While  these  values  are  representative  theycannot 
be  associated  with  a  particular  velocity,  body  shape,  or  vehicle  attitude  because 
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•>  (ktt) 


Figure  8. a  REFRACTION  EFFECTS  OF  FLOW  FIELD  AND  ATMOSPHERf 
LINEAR  DISPLACEMENT  VS  ALTITUDE 
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of  the  involved  interdependence  between  these  quantities.  The  deflection  by 
the  mean  boundary  layer  gradient  is  represented  by  a  curve  based  on  hyper¬ 
sonic  velocities.  However  the  variation  with  velocity  is  not  linear  so  that 
this  curve  is  representative  of  most  practical  vehicles. 
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GLOSSARY  OF  SYMBOLS 


Non-Stereo  Syatema 


NP 

s 

target  offaet  meaaured  along  spherical  earth 

NP' 

t 

target  offset  measured  along  flat  datum  (for  camera) 

NP‘ 

it  _ 

target  offset  measured  along  flat  datum  (for  radar) 

= 

target  offset  error  for  camera 

= 

target  offaet  error  for  radar 

3 

earth  radius 

N 

= 

altitude  of  sensor  above  earth  datum 

& 

= 

angular  offset  of  target 

h 

3 

target  altitude  above  earth  reference  datum 

Y 

= 

target  offset  along  earth  reference  datum 

Stereo  Syatema 

P 

s 

vector  from  base  to  target 

B 

vector  from  base  to  sensor  (  n  ) 

B 

vector  from  sensor  (  n  )  to  target 

T 

(no  aubacrlpt)  =  vector  from  aenaor  1  to  aenaor  2 

= 

the  j  th  measured  variable  m  a  list  of  measured  quantities 

df*) 

S 

the  error  in  the  quantity  called  * 

dx 

s 

combined  film  measurement  errors  (microns) 

dF 

= 

focal  length  measurement  error  (microns) 

dN 

= 

combined  navigation  errors  (all  components  of  error  vectors) 

dSTfS^ 

combined  angular  errors  in  all  orientation  measurement 
(stabilization) 

d/s 

= 

radar  range  measurement  error 

^(h 

= 

error  contribution  due  to  tth  error  source 

^(n) 

a 

altitude  of  sensor  ( n)  above  reference  coordinate  origin 

A 

a 

angular  offset  of  target  from  sensor  1 

^(n) 

a 

focal  length  of  camera  <f>) 

= 

partial  derivation  of  target  position  vector  with  respect  to 
variable  (*). 
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GLOSSARY  OF  SYMBOLS  (Contd.  ) 


c.e.  = 
4 « 

= 

(4) 

LO 

LA 

Ai 

U 

A'.e.c,  = 

a 

f'.o-.E  = 
^  = 

r 

<r* 

r  = 

"X  »> 

5  = 


camera  error 
film  coordinates 

transformation  from  system  j  to  system  ( 
yaw  I  pitch  and  roll  angles 
unit  vector 

change  in  longitude  from  sensor  1  to  sensor  2 
change  in  latitude  from  sensor  1  to  sensor  2 
change  in  flight  heading  from  sensor  1  to  sensor  2 

/E,  -f-e,  ‘T) 

intermediate  quantities  defined  in  Appendix  C  for 
convenience  during  radar  analysis 

film  coordinates  and  focal  length  in  frame  camera 

parameter  for  length  of  camera  to  target  vector 

spherical  coordinates  of  target  with  respect  to  radar  sensor 

residuals  of  least  squares  adjusted  measurements 

vector  of  parameters  for  least  squares  adjustment 

verier  of  parameter,  residuals  after  adjustment 

variance  matrix 

weight  of  the  2  th  measurement 

condition  equation  vector  of  measured  values 

matrices  of  derivatives  of  ^ 

vector  of  Lagrange  multipliers 

weighted  sum  of  squares  of  residuals 


Atmospheric  Effects 

lfi‘  a  tangent  vector  to  light  ray 

n  =  index  of  refraction 

K  a  specific  refractivity 

p  =  density  of  atmosphere 
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GLOSSARY  OF  SYMBOLS  (Contd.  ) 


6  =  angular  displacement  of  ray 

OL  -  ground  displacement  of  ray 

C^c.  =  angle  away  from  vertical 

T  =  temperature 

L  -  lapse  rate  (cLT/dh) 

=  angle  of  incidence  to  boundary  layer 

S  s  deflection  of  ray  vector 

D  -  deflection  modulus 

4  s  argument  of  correlation  function 

^  =  correlation  function 

=  turbulence  scale 

3  s  thickness  of  turbulent  layer 

f'  a  image  displacement 

f  s  focal  length  of  camera 

S(.^  s  distance  (subscript  i'  for  image,  o  for  object) 
A  3  aperture  radius  of  lens 

Ji  =  shimmy  factor 

A  =  altitude 
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APPENDIX  A 


STEREO  PANORAMIC  CAMERA  ERROR  ANALYSIS 

The  triangulation  equations  for  basic  target  location  will  be  derived 
in  this  appendix  for  location  by  two  panoramic  cameras.  Two  frame  cameras 
can  be  handled  by  an  almost  identical  analysis.  The  errors  in  the  target  po¬ 
sition  vector  will  be  evaluated  as  functions  of  the  measurement  errors  by 
differentiation  of  the  target  location  equations. 

A.  1  Derivation  of  Target  Location  Equations 

Each  camera  station  from  which  a  photograph  is  taken  provides 
measurements  of  two  film  coordinates  and  a  focal  length,  which  combine  to 
determine  a  direction  in  space  from  the  camera  to  the  target.  Effectively, 
this  fixes  azimuth  and  zenith  angles  but  without  range  information.  With  two 
photographs,  from  different  stations,  the  target  location  is  over -determined, 
since  four  quantities  are  known  where  three  are  needed.  The  solution  given 
here  for  the  basic  target  location  is  derived  using  equation  4-2  of  Section  4 
(repeated  below)  as  a  starting  point.  The  over  determination  is  eliminated 
mathematically  by  using  only  the  minimum  number  of  equations  sufficient  to 
determine  the  target  location.  Physically,  the  solution  of  these  equations  is 
the  distance  from  one  camera  to  the  point  of  closest  approach  of  the  target 
position  vectors  from  the  two  cameras.  This  point  of  closest  approach  Is, 
of  course,  the  intersection  of  the  lines,  If  they  meet.  Repeating  equation  4-2; 

S  (4-2) 

The  vector  equation,  4-2,  implies  three  scalar  equations.  Rather  than  using 

the  Cartesian  components,  a  better  solution  can  be  obtained  by  taking  scalar 

products  with  the  unit  vectors  determined  by  the  cameras.  Thus,  if  is 

A 

a  unit  vector  whose  direction  Is  fixed  by  camera  number  1,  and  Ag  is  the 
same  for  camera  2,  scalar  products  yield; 

A-1 

^2  •  -4)  «  ^  5  • 
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Since  only  two  of  the  three  possible  equations  are  used,  the  overdetermination 
is  eliminated.  The  magnitude  of  R  may  be  solved  ^rom  equation  A-  1,  as: 


' - - 

The  target  location,  (f?, )  ,  Is  specified  by  the  magnitude,  /?^  ,  and  the 

direction  Rf  with  respect  to  camera  station  1,  as: 

R,  s  R,  R,  A-3 


A.  2  Evaluation  of  Position  Vector  Errors 

The  error  in  Ri  due  to  an  error  in  variable  is  found  by 
differentiating.  The  result  is: 


The  errors  in  the  target  position  vector  can  therefore  be  found 
if  the  errors  in  the  unit  vector  and  the  errors  in  the  vector  magnitude  are 
evaluated. 

A.  3  Errors  in  the  Unit  Vector 

A 

The  unit  vector  Rf  can  be  written  in  terms  of  measured  values 
of  film  coordinates  (  4  ^  )  and  focal  length  {  ■f)  for  the  panoramic  camera 

by  simply  examining  the  geometry  of  the  system,  as  shown  in  figure  A-1. 
Then,  in  terms  of  sensor  1  coordinates: 

«  / 


I 

FatnCa/F) 
F ao9(s/F)\ 
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A 

But  must  be  known  in  earth  fixed  (base)  coordinates.  While  the  base 
coordinates  can  be  chosen  so  that  the  mean  value  of  roll,  pitch  and  yaw  angles 
is  zero,  there  is  error  possible  in  these  angles.  The  transformation  matrix 
must  be  written  out  and  differentiated.  Writing  roll,  pitch  and  yaw  as  p, 
and  the  matrix  is: 


cos^  O  a, 

10  o' 

To,  - 

smp  Cros^  0 

O  1  o 

0  CAe^  -  p 

o  o 

O 

0  »,np  cMp 
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When  only  the  measured  value  is  needed,  it  can  be  taken  as: 


To, 


/  O 
o  t 

o  O 


o 

o 

/ 
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The  error  in  K,  in  base  coordinates  is  then  the  error  in  (To,,R,)  where  R,  is 
measured  in  sensor  1  coordinates.  The  error  in  is  therefore  given  by: 


fjA\  .  T  oL) 

VdA/j/base  /sensor  1 

coord.  coord. 

where  the  products  are  matrix  products. 


A-8 


For  the  error  analysis,  the  measured  variables  for  the  system 


must  be  listed: 

'  li  1 

•  Pi 

1 

/f/j,  » 

>  sensor 

1  orientation 

^  J 

•  Sk  St» 

1 

'  Sy  Siy 

^  sensor 

2  orientation 

^  S|B 

I 

/t/fl, «.  /* 

] 

rt/g  »  Siy 

>  sensor 

2 

/j/j  »  Sig 

1 
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The  variables  -t/}..#-. 4  are  dependent  on  how  the  navigation  veii-or  v.  measured. 
If  both  sensors  are  mounted  on  the  same  aerial  platform,  and  the  two  sensor 
records  are  made  during  the  same  flight,  then  the  navigation  vector,  S 
involves  errors  due  only  to  the  change  of  position  during  the  interval  between 
records.  In  that  case,  components  of  S  If  the  sensor 

records  are  made  from  different  platforms  and  therefore  different  flights, 
the  vector  s  involves  errors  from  the  two  navigation  vectors  S^.  and  S, 

In  that  case  are  the  components  of  Sj  ,  For  one  flight  5  is 

independent  of  errors  in  Sf  ,  while  for  two  flights,  S  is  given  by; 

(4-3) 


It  is  convenient  to  differentiate  always  with  respect  to  components  of  "S' 
and  take  the  change  of  dependent  variables  into  account  separately  as: 


i  (function]  d  (function) 
15/  ■  35/ 
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^  (function)  ^  (function) 


aCSt)/ 

9  (function) 


dfS,)/ 


<95,' 

3  (function) 
35/ 


AAi 


d(5,)/ 
35 


for  2  flights 


where 


35/ 


s  1,  and 


^S, 


i(S, 

I 

=  -1  as  follows,  since 


S  »  5,-  S. 
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Therefore,  derivatives  with  respect  to  S/  will  be  written  as  if  ^.r.k 
were  the  components  5/  and  derivatives  with  respect  to  will  be 

omitted  from  the  list.  The  values  of  unit  vector  error  may  be  computed  in 
base  coordinates  by  differentiating  equations  A- 5  and  A- 6  and  combining  accord 
ing  to  equation  A-8.  Omitting  zero  values,  these  are  ( ( ,  a  ,  and  Prefer  to 
sensor  1): 


r 


I 

(/*•*■  F*)'"* 


A- 11 
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^(R,)x 

FI 

9(iU  ^ 

-  P  coi(a/F) 

^'^11 

(1^4- F*)'^' 

aC.^/.Vit 

+  ^  (sIf) 

^n/if 

di^,)y  i 

IF  am  (a/p) 

d(R,)y 

-COS  Cs/rJ 

at^p»yit  * 

dn\ 

* 

Mix  \ 

• 

/  y#-»\  /  Vj»  \  ^  /j»\  *_ 

r‘am(s/F)  j 

•  / 

aA^,  ( 

(9/F)ca»  (s/p)  -  at'>(s/F) 

a'4/«'*)  j 

Fco»  (alF) 

a(«,)y 

-/ 

a  AT,, 

JPTp^  » 

d^,t 

Mix 

IF ooi  Ca/p) 

Mi,)» 

tin  (t/F) 

a-vj 

UUF')^ 

Mix 

d/\^t 

~{cx>9(9/F)  *-Ca/F)9m(9/F)]+ 

Mix  . 

1 

u 

af/ej,  ^ 

1 

a^cr. 

dA/}, 

(TsT^tyv, 
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A.  4  Errors  in  Position  Vector  Magnitude 


The  values  of  [  hRf  !  d't/'j-  ]  are  derived  by  differentiating  the 

equation  for  Rf  (A-2)  with  respect  to  each  Equation  A- 2  is  rewritten 

with  quantities  grouped  for  convenience; 


where: 


R,  e 


A-BE 


A  = 

A  A 

f » 


Then,  the  derivatives  of  R/  are: 


Where, 

dA 


__i _ r  hA  _  dB 


(zee -B) 


I 


bb  _  ^  a^i  ^ 


A-12 


A-13 


A-14 


A-15 


A-16 


With  the  above  grouping,  A  ,  3  ,  and  E  and  their  derivatives  are 

scalars.  However,  the  individual  terms  in  the  last  three  equations  contain 
scalar  products  of  vectors.  Therefore,  all  of  these  vectors  must  be  trans¬ 
formed  to  the  same  coordinate  system.  The  logical  system  to  choose  is  an 
earth  fixed  system.  Since  the  target  is  being  located  by  evaluating  Rf  , 
the  earth  fixed  system  is  chosen  at  sensor  1.  These  coordinates  are  the 
translated  base  coordinates  already  implicitly  used  in  section 4,  as  well  as 
for  the  computation  of  [  dR^  / ].  Thus,  R,  has  been  transformed 


to  this  coordinate  system,  and  S  is  measured  in  this  system  (any  earth 
fixed  system  would  do).  Only  must  be  transformed.  This  transforma¬ 

tion  is  done  in  two  steps: 

1.  Transform  from  sensor  2  coordinates  to  earth  fixed  coordinates 

at  sensor  2.  This  transformation  is  exactly  to  form  of  that  used 
for  ,  and  the  values  of  [  ]  can  be  used  for 

if  all  reference  to  sensor  1  is  replaced  by  sensor  2.  Therefore, 
a  and  -f  are  used  from  sensor  2,  and  ^  to  5,  /Oto  n) 
are  replaced  by  to  fs) 

2.  Transform  from  earth  fixed  coordinates  at  sensor  2  to  earth 
fixed  coordinates  at  sensor  1.  (The  errors  involved  are  due  to 
navigation. ) 

The  second  transformation  matrix  depends  on  the  relative  orien¬ 
tation  of  the  earth  fixed  systems  at  sensor  1  and  sensor  2.  The  different 
orientation  is  due  to  the  possibility  of  different  flight  directions,  and  to  the 
change  of  direction  of  vertical  across  the  surface  of  the  earth.  If  the  X 
direction  at  sensor  1  is  called  north,  then  a  latitude  and  longitude  can  be 
connected  to  the  components  of  the  navigation  vector  using  the  usual  defini¬ 
tions.  The  change  in  flight  direction  will  be  measured  as  an  azimuth  angle. 
These  angles  are; 


longitude  a  LO  *  Sv/^« 

latitude  a  LA  *  A- 17 


azimuth  a  A£  a  flight  direction  of  sensor  2  with  respect 
to  1)  where  Sg  is  the  radius  of  the  earth.  Using  {  C  )  to  mean  (cos)  and 
(5  )  to  mean  (sin),  the  transformation  matrix  is: 


~(SM)(aOh(CAi)(SLA)(SLO) 
•(^KlXlSLO)  -  (CrnSLAYp-O) 


(sM)(af<) 

(CMXCLO)  *  (SM)(SLA)lbLO) 
(CAi)(iL0)~(SMX9LA)(CL0) 


(SLA) 

-COaKslo) 

(Ci.A)(CLO) 
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Since  error  in  flight  direction  has  been  taken  into  account  in  the  first  trans¬ 
formation,  the  value  of  AB  is  known  in  this  transformation,  without  error. 

A 

The  error  in  Rt  is  found  in  translated  base  coordinates  as: 


A-19 


The  terms  in  this  equation  have  been  computed,  except  for  ( dT^g/d^^^  . 
The  terms  can  be  found: 


Tat  in  equation  A- 18 

use  equation  A- 1 1,  changing  to  Rg  ,  and 

changing  to  j  •  n, 

^  J  »  /O,  ft,  tZ  to  j  /J,  /4^  IS 

use  equation  A- 5,  changing  k,  to  I  to  li.SgJi 
The  term  [  ^  Tos  ]  contains  errors  due  only  to  Sx  &nd  Sy  (variables 
ntg  and  ).  These  are  (using  the  same  notation  as  equation  A- 18): 


f-i'OiiX'SU)  -fSAEXstA)  (CLA) 

(UaKCLAySLO)  (SAi)(CLA)('iLO)  <5LA)(tL0) 

L-^CA2;(CtAKCl-0)  -(5A2)(dtA)(CLO)  -(SLAKclo) 
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d  Tog 
d/t/j. 


O  O 

-(SAi)(SLO)  -(CAl)(SLA)(CLO)  (CAi)(Sl0)  -(SAlXSLf^CLO) 
[  (SAlXCLO)  -  fCA»)fSLA)(SL6)  -(CAlKOO) -(bAlXSiA)(CLO) 


O 

(CLAXCtO) 

(CU)(^LO) 
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Because  of  the  matrix  multiplications  for  coordinate  Ira  iisfor <naticn,  ths 
equations  would  be  overly  lengthy  if  expressions  for  [dA,  / dirj  ]  were  written 
out  explicitly.  In  practice,  when  these  values  are  computed,  they  are  evalu¬ 
ated  step  by  step,  using  the  expressions  already  given.  Equation  A- 19,  with 
the  above  expressions,  provides  a  value  of[  ]  and  Rg  . 
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A  ^ 

[  found  from  Equations  A-ll  and  A-5,  [  ] 

can  be  evaluated  using  Equations  A- 13  to  A- 16,  plus: 


S  ■  [Sk  .  Sy  ,  5«] 


and  (omitting  zero  values): 


i£ 


« 


[  1.0.0] 


as 

* 


[o./.o] 


[0.0.  i] 
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A.  5  Results 

A  4 

The  values  of  [  ],  ,  and  [  /dvj  ] 

have  been  derived  for  photographic  stereo  in  this  appendix.  These  terms  are 
used  to  provide  the  required  list  of  values  of  [  ]  for  Section  4. 
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APPENDIX  B 


PANORAMIC  CAMERA  WITH  RADAR,  ERROR  ANALYSIS 


In  this  appendix,  the  target  location  equations  are  derived  for  a  panoramic 
camera  combined  with  a  side-looking  radar.  The  values  of  [  3(Ri)  / dfifj  ] 
are  derived  by  differentiation  of  the  location  equations. 


B.  1  Derivation  of  the  Location  Equations 

The  target  is  located  by  a  combination  of  directional  information 
provided  by  the  camera,  and  range  information  provided  by  the  radar.  Since 
the  two  sensors  are  not  necessarily  at  the  same  point,  the  solution  is  not  the 
trivial,  separated  form  given  by  that  special  case  I  oince  aftiii'iuth  information 
provided  by  radar  is  relatively  low  precision,  only  the  range  measurement  is 
used  in  the  basic  location  analysis.  The  solution  therefore  takes  the  form  of 
finding  the  intersections  of  a  sphere  and  a  line,  the  sphere  being  the  surface 
of  constant  range  from  the  radar  and  the  line  being  the  direction  fixed  by  the 
camera.  Equation  4-2  of  Section  4  is  the  basis  for  the  derlvatfon  of  the  location 
equations.  Repeating  Equation  4-2; 


R, 


5  +  R, 


(4-2) 


The  direction  of  Ri  is  known  from  camera  measurements  as  given  in  Appendix 
A.  Repeating  the  value  of  this  unit  vector  (R|)  for  panoramic  camera: 


where:  s  , 


are  film  coordinates 


I 

F  5m 
,FMn*/F  . 


(A-5) 


F  a  focal  length 

In  Equation  4-2,  repeated  above,  the  direction  of  R,  is  known  5/  , 

the  navigation  vector,  is  known;  the  magnitude  of  is  measured  by  radar; 

and  the  magnitude  of  ^  is  required.  The  magnitude  of  R,  may  be  found  by 
taking  the  scalar  product  of  R,  with  each  term  of  4-2: 


R‘  e  R,  •  +  R,  ‘R, 
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B-2 


Similarly,  the  scalar  product  of  Equation  4-2  with  3  and  Rg  yield: 

??r**5*“ =  S‘  +  ^*3" 

B-3 


Combining : 

=  2(kr'T)R,-h(R*-5*)  B-4 

Equation  B-4  is  a  quadratic  scalar  equation  for  the  magnitude  of  Rf  . 
Solving  for  Rf  : 

S,  a  J  f  -/(R  ‘5  -3*)  B-5 

The  two  solutions  correspond  to  the  two  intersections  of  tlie  line  with  the 
sphere.  This  ambiguity  will  be  treated  leie  r  All  uie  neCvion*  T'aao  Qu&ntity 
^  can  be  evaluated: 

R,S^(RX  S*+('R,)y  Sy*(RX^^ 

•  ~yp7f^j^ +  cos  i/pj 
It  is  also  convenient  to  define: 


U  =  1 'T)  *  •<•  (  R,*- S*  ) 

B-7 

so  that: 

II 

c 

B-8 

B.  2  Error  Analysis. 

-ix  A 

The  vector  giving  target  location  is  R,  specified  as  R,  R, 

by  Equations  A-5  and  B-8.  Except  that  R,  is  computed  by  a  different  equation 

for  the  panoramic  camera-radar  combination  from  that  for  two  cameras,  the 

A 

analysis  is  identical  to  that  given  in  Appendix  A.  The  values  of  Rf  and 
can  be  taken  from  Appendix  A.  Only  the  [BRi/dnfjJ  must  be 
recalculated  and  a  slightly  revised  list  of  variables  is  needed  to  account  for 
sensor  2  now  being  a  radar  unit.  The  new  variable  list  is: 

,y  s  1  to  12  are  the  same  as  in  Appendix  A. 

,  the  radar  range  measurement. 
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The  same  c or.aide ration  as  described  in  Appendix  A  must  be  used  when  navigation 
measurements  are  treated,  with  respect  to  single  or  separate  flights  over  the 
target  area. 


The  foregoing  list  of  quantities  available  in  Appendix  A,  plus  the 
values  of  are  needed  for  Equation  B-9.  In  order  to  differentiate 

it  is  convenient  to  evaluate  5)/d/»/j]  and  separately. 

They  are: 


dU 

dnfl 


d(K..T) 

'd's  \ 

B-9 

,/ 

/ 

d(R|‘S)  dRz 

.  ] 

B-10 

”  U 

In  B-9.  the  first  term  is  listed  in  Appendix  Ai  only  the  scalar  product  needs 
to  be  carried  out.  The  second  term  is  only  non-zero  for 
The  values  of  ]  can  be  tabulated; 


dR 


H 


i(r,'T)  ^ 


Ss  iinCi/F)  "  5y  coi( i/F) 


B-11 


d/e, 


ijMl 


A  ♦ 


ag/  jFsm(*/F)  I 

o**F*y^  U 


P3>nG/F) 


_ y  Fcoa(»/F)  .  I  (R,-^)  Pco3(i/F)  1( 
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see  Appendix  A 


ifZ,  ^  ^  dB.,  ] 

dnJ-^  dn/'y  j 

3/g«  f  Sy  Paos  ^ s/F)  -  S,  F3//>  fs/r)]  [  I  ^  /^|  *5  j 

1 - - IL  i  u  / 

«3«.  _  [5.2  -SKFco4(9/F)ir.  /ie, 'Sill 

a«-«  L — L'  (“0“/] 

rSvF9in^e/F)-5ynr 

a^s  [  r/**F‘F-  lr(  a  jj 

a>w;j 


Since  would  be  given,  in  an  error  analysis,  and  S  !> 
known,  the  value  of  U  may  be  found  from: 


B-12 

thereby  eliminating  the  ambiguity  originally  in  the  sign  of  U  ,  The  sets  of 
derivatives  given  by  equations  here  and  in  Appendix  A  yield  the  desired  results 
for  [d^,/avyj  . 
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APPENDIX  C 


RADAR  STEREO  ERROR  ANALYSIS 


In  Section  4,  the  format  of  the  error  analysis  is  laid  out.  The  purpose 
of  this  appendix  is  to  establish  the  form  of  the  basic  target  location  equations 
for  radar  stereo  and  to  evaluate  f  j  for  this  case. 

C.  1  Location  Equations 

Each  radar  sensor  is  assumed  to  be  a  side-looking  system  with 
no  elevation  measurement.  The  range  measurement  is  usually  the  most  accurate 
measurement,  while  the  azimuthal  position  is  within  the  spread  of  the  antenna 
pattern,  which  is  assumed  here  to  be  centered  at  a  right  angle  to  the  flight 
line.  Two  measurements  (a  range  and  azimuth  determination)  from  each  radar 
sensor  overdetermine s  the  target  position.  A  range  from  each  sensor  plus  one 
azimuth  angle  will  locate  the  target,  then  using  the  other  azimuth  angle  would 
reduce  the  location  error  (azimuth  error  would  be  reduced  approximately  by 
Vi  )•  Therefore,  the  basic  system  will  be  chosen  to  use  two  ranges  and  one 
azimuth  angle,  and  the  values  assigned  to  azimuth  error  will  be  assumed  to 
take  account  of  the  factor  of  YF  . 

The  required  location  equations  are  the  equations  of  2  spheres  and 
a  plane  since  two  ranges  are  measured  and  the  antenna  beam  is  essentially 
planar;  the  intersection  point  being  the  target  location.  Referring  to  Figure 
4.  1,  Section  4,  and  writing  the  equations  in  sensor  coordinates: 

{K,)*  s  X* 

X,  •  0 

where  subscript  1  means  "measured  in  coordinates  of  sensor  1.  " 

Rotating  the  coordinate  system  such  that  the  new  X  axis  is  in  the 
direction  of  the  base  X  axis  provides  results  directly  in  base  coordinates, 
except  for  translation.  (It  is  again  convenient  to  assume  the  base  coordinates 
are  oriented  with  £  being  vertical  at  sensor  1,  thus  allowing  earth  curvature 
to  be  ignored  for  the  error  analysis.)  The  equations  for  spheres  remain  the 
same,  except  that  quantities  are  now  measured  in  translated  base  coordinates. 
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The  transformation  for  X,  i>: 


X,  -  XCM  tp  -  Vs/rt  ^ 


C-2 


where  p  s  angle  between  flight  direction  at  sensor  station  1  and  x  axis  of 
base  system.  The  new  equations  are  therefore: 

R*‘  •  Y* 

•^x-sj‘+ry-Sy)*+^#-5,)*  c-3 

A'x-fly  a  o 

where  ,  These  three  equations  can  be  solved  for 

by  simple  algebra,  giving: 


where: 


«  M.Ss  ±  CD 

*•  — e - 


AX  » 


C  a 


^/‘Ti?aV3* 

2 

5-5*  ^  A'Sy 


C-4 

C-5 

C-6 


C-7 


The  plus  and  minus  choice  in  the  value  of  B  is  cbe  to  the  two  possible  inter¬ 
sections  of  the  plane  and  two  spheres.  In  the  error  analysis  this  will  be  re¬ 
solved  because  target  location  is  an  input,  in  practice  the  two  2  values 
would  be  sufficiently  different  as  to  be  separable  by  approximate  altlti^do 
measurement.  (One  intersection  may  well  be  above  the  sensors.) 

For  the  error  analysis,  XiY,  2  must  be  differentiated  with  respect 
to  each  measured  variable.  In  this  case,  the  variables  are: 


•  Sf  Si« 

^  •  Eg 

^  •  S,a 

/*5  -  f 

^  »  S|y 

C-8 

*  Sx  Si* 

'  S)» 

5/  ^ay 


The  navigation  variables  ^  are  treated  Just  as  in 

Appendices  A  and  B,  with  respect  to  one  or  two  flights  over  the  target.  There¬ 
fore,  only  derivatives  with  respect  to  variables  1  to  6  are  required. 
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C.  2.  Error  Analysis 


Differentiating  Equations  C-4  and  C-5  for  X  and  Y  : 

dir^ 


c^/i^  d\/j  diJ2 


C-10 


These  expressions  will  be  evaluated  for  particular  j  values  after 
the  values  are  derived. 

The  derivatives  of  A,  fi  ,  C  ^  D  ,  £  ,  and  hA  are  necessary  to  evaluate 
.  These  derivatives  are  listed  below,  with  zero  values  omitted. 


alr^  "  *  ^ 

diM  _  dL  M  _  clM  _ 

17, 

dC  cLC  _  oLC  , 

.S  ,  ^  ■  A 

dv,  '  dir.  ‘77  ‘17 


dD  .  et(e-M) 


dp  _  MSm 
^  di/'i  “  P 


dO  ^  g*C  ^  ^  ^  ) 

dv-j  O  c^vj  *  P 

,  3s  rg*“  A/; 

c^v;.  ^  ^  p 


C-ll 


In  terms  of  the  above: 


£  j  dvj  *  dvj  dv^  “  di/j  ^  duy  j 
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After  all  differentiation  has  been  completed,  no  generality  is  lost 
in  the  error  study  if  ip  la  taken  as  -  0.  This  can  be  done  because  no  direction 
has  been  set  for  the  axis  of  the  base  coordinates.  With  this  simplification: 

a'--  /  ,  &  ‘O  C-13 


The  values  of  -j —  . -  ,  and  -5 —  may  be  tabulated.  They 

'd-O'J  dvj  ^ 

are  listed  below,  with  zero  values  omitted.  Values  are  found  by  substitution 
of  C-11  and  C-13  into  Equations  C-9>  C-10,  and  C*  12. 


dx 


'  r 


C-14 


du;  Y  Y  dv;  *  du*  Y  dix/ 


J  •  2  to  (t 


di  R, 


Sy(£‘Mjl 

-IT- I 


^  El! 

dtx,’  F  1 


S.t 


- 

dt/ir 


Throughout  Equations  C-14,  the  double  valued  character  of  -i 
causes  doubled  values  in  dt/dvj  ,  This  ambiguity  can  be  removed  in  the 
error  analysis,  since  f  is  an  input  parameter.  If  it  is  noted  that  p  appears 
with  every  plus-minus  pair,  the  source  of  the  ambiguity  may  be  found,  since 
t>  was  solved  for  as  a  square  root  and  arbitrarily  called  positive.  For  the 
error  analysis,  0  can  be  defined  equivalently  by  C-6,  as: 


The  expression  tp  may  be  replaced  by  the  equivalent  form  which  will  assign 
the  proper  sign  to  the  quantity,  without  ambiguity.  Substitution  of  the  values 
ol  C  ,  E  and  M ,  and  algebra  will  simplify  C-15  to: 
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C-16 


Using  Equation  C-16.  the  solutions  for  dii/cLv^  are  unique. 

Since 

i  C-17 

the  values  of  d(R^ )/ cLy]  are  listed  as  Equation  C-14,  completing  the 

derivation. 
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APPENDIX  D 


EFFECTS  OF  MULTIPLE  COVERAGE  OF  A  SINGLE  TARGET 


In  many  cases,  several  ^-^’'nsor  records  c^'iitain  data  on  the  location  of  a 
target.  Since  the  target  can  be  located  with  respect  to  one  sorsor  by  meas¬ 
urement  of  three  independent  quantities,  it  is  both  possible  and  often  probable 
that  excess  or  redundant  measurements  will  be  available.  Redundant  meas¬ 
urements  can  be  used  to  increase  target  location  accuracy  by  adjusting  the 
measured  quantities  according  to  a  criterion  such  as  maximum  likelihood. 

Weighted  least  squaves  analyses  are  available  in  the  literature  which  are 
concerned  with  adjustment  of  measurements  for  photogrammetrlc  control 
extension.  D.  C.  Brawn's,  Reference  8,  report  contains  a  relatively  complete 
treatment  of  the  theory  of  weighted  least  squares.  The  method  presented 
there  will  be  used  as  a  basis  for  the  treatment  of  adjustments  given  here. 

The  form  of  the  least  squares  analysis  can  be  forced  to  remain  fixed, 
independent  of  what  sensor  combinations  are  used  to  locate  the  target,  and 
how  many  sensor  records  contain  information  about  the  target.  The  problem 
can  be  reduced  to  one  of  writing  the  condition  equations  in  proper  form.  In 
all  cases,  the  condition  equations  will  be  written  (using  the  same  notation  as 
in  previous  sections)  as: 

W,  •  Rt  *■  5  D  1 

where  and  Rt  are  vectors  from  two  sensor  stations  to  the  target,  and  S 
is  the  navigation  vector  giving  the  location  of  sensor  2  with  respect  to  sensor 
1.  Of  course,  all  vectors  must  be  specified  in  a  common  ''oordlnate  system, 
so  that  some  vectors  may  first  have  to  be  transformed  tu  the  proper  system 
before  being  used. 

The  firt^t  pair  of  stations  provides  one  vector  condition  equation  (e.  g.  , 

3  scalar  equations).  Each  additional  station  provides  another  vector  equation, 
since  each  such  equation  expresses  the  fact  that  the  vector  from  one  station 
intersects  the  vector  from  station  1,  at  the  target. 
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In  order  to  carry  out  the  analysis  with  uncorrelated  measurements,  the 
form  in  which  each  vector  is  expressed  depends  on  what  sensor  is  considered, 
since  the  vector  must  be  expressed  in  measured  quantities.  For  example, 
the  vector  for  a  frame  camera  would  be: 


where; 


f 

CL 


R,  «  a 


are  film  coordinates 
is  focal  length 

is  a  multiplier  which  is  unknown,  since  the  camera 
does  not  measure  distance  to  a  target 


D-2 


For  a  side-looking  radar  system,  however,  the  vector  would  be: 


where; 


»  r 


t6%  c* 
art  r  LOh  S 

S4I^9“  JUm  P 


r>-3 


r  is  the  measured  range 

cr  is  the  antenna  squint  angle  (normally  W/t  for  side-looking) 
e  is  the  elevation  antle  (usually  not  measured) 


These  vectors  are  expressed  in  sensor  coordinates,  and  would  have  to  be 
transformed  to  earth  coordinates  by  a  matrix  involving  roll,  pitch  and  yaw 
angles  of  the  sensor. 


The  values  of  the  quantities  used  would  be  measured  values  for  all 
measured  quantities,  and  estimates  for  unmeasured  quantities.  The  estimates, 
would  be  obtained  from  the  basic  location  equations  derived  in  the  previous 
section*  • 


The  weighted  least  squares  analysis  can  be  carried  out  by  following  the 
form  given  in  Appendix  A  of  Reference  8.  Since  the  theoretical  aspects  of 
the  analysis  are  well  covered  in  the  reference,  only  a  partial  treatment  is 
given  here,  in  order  to  indicate  the  use  of  these  procedures  with  sensors  other 
than  cameras.  A  brief,  simplified  example  is  given.  The  simplest  possible 
case  of  the  use  of  redundant  data  will  be  treated,  in  order  to  show  how  sensors 
other  than  cameras  can  be  used  for  target  location. 
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The  case  to  be  described  is  tarjjet  location  by  means  of  one  vertical  frame 
camera  at  one  station^  and  a  aide*looking  radar  at  a  nearby  station.  Since  the 
camera  determines  the  direction  of  a  vector  and  the  radar  determines  a  magni¬ 
tude  and  one  direction  of  a  second  vector,  four  quantities  are  known,  so  the 
target  is  overdetermined  by  one  measurement.  Navigation  and  orientation 
measurements  could  be  included  in  the  adjustment,  but  this  only  lengthens  the 
equations  without  adding  to  the  illustrative  value  of  the  example.  Moreover, 
the  orientation  (roll,  pitch,  and  yaw)  angles  can  be  set  as  zero,  thus  making 
all  transformation  matrices  simply  unit  matrices.  Wherever  possible,  the 
notation  of  Reference  8  will  be  used,  so  that  comparisons  can  be  made  easily. 


It  will  be  assumed  that  an  approximate  target  location  has  been  found, 
so  that  rough  values  are  available  for  all  quantities.  The  first  problem,  then, 
is  to  write  the  condition  equations.  These  are  given  above.  The  form  of 
and  ffg  is  taken  from  Equation  D-2  and  Equation  D>3.  The  condition  equations, 
in  vector  form,  are: 


»  < 

s  < 

.  * 

ees  r 

Sa 

C6 

Yw 

-f 

i<n«r  sesfi 

- 

Sy 

4 

4 

er  E 

1  m 

Si 

In  this  equatlen,  the  components  of  S  are  assumed  to  be  known,  not  to  be 
adjusted,  and  therefore  constant.  Because  of  the  assumptions  of  known  orienta¬ 
tion,  no  transformation  matrices  appear.  The  quantities  to  be  adjusted  can  be 
separated  into  2  groups:  those  measured  quantities  (  Xj)  whose  weighted 
adjustments  are  to  be  minimized,  and  the  unmeasured  quantities  (/^)  which 
are  adjusted  with  zero  weight  (called  parameters),  These  are: 


x-f 

f 

r 

e- 
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Labelling  measured  or  estimated  values  with  a  superscript  zero,  the  vectors 
of  correction  values  may  be  defined: 


^  -It* 

7* '  ?  - 


D-6 
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The  5  measured  quantities  are  assigned  weights  which  are  inversely  pro¬ 
portional  to  the  variances  of  the  measurements,  and  a  variance  matrix  is 
written  with  the  variance  of  the  J  th  measurement  as  ihejj  th  ms^rix 
element: 


O'’  a 


O 

o 

o 
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-/ 
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o 

0 

o 


D-7 


The  condition  equations  {D-4)  are  linearised  by  expanding  in  a  Taylor  aeries 
about  the  measured  values,  keeping  only  linear  terms.  It  is  convenient  to 
define: 


£  =  the  vector  condition  equation  with  measured  values  inserted 

in  place  of  true  values.  (  I  th  component  is  €|'  ) 


“v  ■  -  *“  * 

At  6  °  matrices  with  elements  .  biit 
The  condition  equation  (D-4)  may  be  written  in  linearised  form  as: 

*■  OS  +6*0  D‘9 

For  the  present  case,  the  matrices  A  end  5  are: 


'off* 

0 

0 

0 

r* 

1 

a* 

o 

A  » 

0 

a* 

0 

•coaE* 

0 

B  m 

0 

0 

a* 

0 

f*  -r'aesF* 

« 

Using  the  method  of  Lagrange  multipliers,  the  weighted  sum  of  squares  of 
residuals  in  within: 

5  D-11 

where  the  vector  of  Lagrange  multipliers  is:  X|)  for  this  case. 

The  magnitude  of  S  is  to  be  minimised.  This  process  is  carried  out 
by  differentiating  S  with  respect  to  each  free  variable  (i.  e.  ,  all  &nd 

Sa  ).  Each  partial  derivative  is  set  eqxial  to  zero  to  obtain  a  minimvun 
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of  5  .  Differentiating  D-ll  and  writing  the  results  in  matrix  form  gives: 

D-12 

-  O  d-13 

where  D«12  comes  from  differentiating  with  respect  to  ^  and  D-13  from  6  . 
For  the  example  given,  Equations  D-9,  D-12,  D-13  represent  a  set  of  10 
equations  (3  +  5  +  2)  in  10  unknowns  2  &'s  ,  S  7.’»)  to  be  solved. 

Since  the  purpose  of  this  example  is  to  outline  the  procedure  to  be  used, 
rather  than  to  solve  a  particular  problem,  the  straightforward  but  extremely 
lengthy  solution  will  not  be  completed.  The  methods  of  solution  and  error 
analysis  can  be  found  in  Reference  1.  All  of  the  required  matrices  have 
been  set  up  for  the  problem  considered,  thus  showing  how  such  a  sensor 
combination  as  camera- radar  would  be  treated. 
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